WO 03/081484 




PCT/DE03/01002 



1 

Verfahren und Vorrichtung zur Ermittlung von De format ionen 



Beschreibung: 

5 

Die Erfindung betrifft ein Verfahren zur Ermittlung von 
Deformat ionen, wobei diese erfasst werden und wobei untersucht 
wird, wie eine Zufuhr von Energie, hervorgeruf en durch eine 
Belastung, einen Korper verf ormt . 

10 

Die Erfindung betrifft ferner eine fur die Durchfuhrung des 
Verfahrens geeignete Vorrichtung. 

Es sind verschiedene Verfahren bekannt, mit denen die 
15 Deformation eines Korpers unter der Einwirkung von Kraften 
ermittelt wird. 

Die bekannten Verfahren sehen vor, de:. .KSrper durch ein 
Gittermodell darzustellen und fur jeden einzelnen Gitterpunkt 
20 einzeln eine zu erwartende Anderung der relativen Position bei 
der Einwirkung einer Kraft zu ermitteln. 

Vor allem fur den Fall, dass der Korper eine komplexe Gestalt 
aufweist, sind die bekannten Verfahren sehr aufwendig. 

25 

Beispielsweise erfordert die Ermittlung der Deformation eines 
Fahrzeuges bei einem Verkehrsunf all einen sehr hohen 
Berechnungsauf wand . 

3 0 Obwohl die in . den bekannten Verfahren beinhalteten 
mathematischen Modelle einfach sind - die zugrunde liegenden 
mathematischen Modelle beruhen im Wesentlichen auf Arbeiten 
Eulers aus dem 18. Jahrhundert - sind bisher keine Verfahren 
bekannt, die sich prinzipiell von den auf Euler beruhenden 

35 Verfahren unterscheiden. 
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Der Erfinder hat daruberhinausgehende theoretische Ansatze zur 
Deformationstheorie entwickelt und in einer Reihe frei 
zuganglicher Artikel dargestellt. 

5 

Der erste dieser Artikel „An approach to deformation theory 
based on Boyle's law. I. The concept, and kinematics of pure 
shear" behandelt ein Konzept einer Kinematik von Festkorper- 
deformation auf Grundlage einer Adaption der idealen 
10 Zustandsgleichung pV = nRT. Durch diese Theorie werden sowohl 
isotrope als auch anisotrope Materiale und Beanspruchungen 
erf asst . 

Weiterentwicklungen zu der Kinematik sind in dem Artikel „An 
15 approach to deformation theory based on Boyle's law. II. 
Kinematics of simple she--." and some energetic considerations w 
dargestellt. 

Das Erfordernis der Berucksichtigung einer umfassenden 
2 0 Betrachtung auch geometrisch sehr einfacher Def ormationen in 
grundsatzlich alien 3 Raumdimensionen wurde in dem Artikel 
dargestellt: „An approach to deformation theory based on 
Boyle's law. III. Three-dimensional properties of plane-strain 
deformation: origin of conjugate joint sets, sheath folds in 

2 5 plastic shear zones, and turbulence in viscous flow". 

Der Erfindung liegt die Aufgabe zugrunde, ein gattungsgemaSes 
Verfahren so weiter zu entwickeln, dass sich mit ihm 
Verformungen moglichst schnell und zuverlassig ermitteln 
30 lassen. 

Erf indungsgemaS wird diese Aufgabe dadurch gelost, dass ein 
jeglicher Punkt Q in einem Festkorper als Ref erenzpunkt gedacht 
werden kann, der in einer Norm-Einheitsentf ernung in alien 

3 5 Richtungen von Korrespondenzpunkten P umgeben ist. Diese 
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Korrespondenzpunkte P dienen als Richtungsreprasentanten aus 
der Sicht von Q und formen in ihrer Gesamtheit eine Oberflache 
von idealiter kugelf ormiger Gestalt. Die Ortsanderung aller P 
relativ zu Q als Funktion der wirkenden Krafte beschreiben die 
5 Deformation am Punkt Q. Im herkommlichen Verfahren wird ein 
Punktegitter vordef iniert , und alle Punkte dienen gleichzeitig 
als Bezugspunkt und Korrespondenzpunkt . Dadurch ist die 
Genauigkeit des Verfahrens abhangig von der Zahl der 
Punkt epaare, die zwar sehr hoch, aber endlich und 

10 notwendigerweise liickenhaft ist. Im Gegensatz dazu kann beim 
neuen Verfahren jeder Punkt als Bezugspunkt Qi dienen, aber 
damit ist die entsprechende Oberflache A aller zugeordneten Pi 
fur dieses Qi festgelegt; dass ein Punkt Pi fur sich als ein 
Punkt Q 2 fungieren kann, ist aus der Sicht von Qi unerheblich. 

15 Auf diese Weise gibt es fur jedes Q eine unendlich groSe Zahl 
von P, '-ur;d innerhalb des Gultigkeitsbereiches des Ansatzes eine 
■ unep^J -':L jroBe Zahl von Q, was eine Darstellung als homogener 
Dif r^rftt^^iatisatz erst moglich macht . 



Bisheriges Verfahren 


Hier vorgestelltes Verfahren 


Das Netz von Knotenpunkten muss 
erst definiert werden (Abb. 01) . 


Jeder Punkt Q ist von 
vornherein mit einem Feld von 
Punkt en P vorhanden (Abb. 01) . 


Die Zahl der Knoten im Netz ist 
zwar beliebig groS, jedoch 
endlich. 


Die Zahl der Punkte P is 
unendlich. 


Die Zahl der Verbindungslinien 
zwischen Punktpaaren wachst mit 
der Zahl der vorhandenen 
Punkte, ist aber im Prinzip 
diskret, erfordert daher 
jeweils eine eigene 
mathematische Modellierung . 


Die Zahl der Verbindungslinien 
zwischen einem Q und alien 
zugeordneten P ist unendlich. 
Ihre mathematische Behandlung 
kann daher kontinuierlich sein. 
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Die Strecke zwischen zwei 
Knoten ist finit. 


Die Strecke Q^P ist 
maSstabunabhangig . 


Zwei beliebige Knoten im Netz 
dienen als Anfangs- und 
Endpunkt einer gegebenen 
Strecke . 


Die Strecke r 0 ist die 
Verbindungsline zwischen Q und 
P, aber nicht die zwischen zwei 
Q. Es besteht keine Beziehung 
zwischen dem Oberf lachenpunkt 
P x eines Qi zugeordneten 
Normvolumens und einem anderen 
Q 2 , (auch dann nicht, wenn Pi 
gleichzeitig ein P 2 sei,) es 
sei denn, ein Pi falle zufallig 
auf ein Q 2 . 


Deformation wird modelliert 
durch Veranderung dieser 
Strecke in Lange und/oder 
Richtung relativ zu dem einen 
oder dem anderen Knoten, und 
dann relativ zu einem externen 
Bezugspunkt . 


Deformation wird modelliert mi t: ? 

l 

Hilfe von r relativ nur '-*u\? 

f • 

seinem eigenen Bezugspunkt Q, 
und erst indirekt mit Bezug zu 
einem externen Bezugspunkt, 


Die Def ormat ionsmodel 1 ierung 
bezieht sich auf die x 
diskreten Verbindungslinien 
aller Punktpaare. 


Die Def ormat ionsmodell ierung 
geschieht homogen- 
kontinuierlich uber alle P 
relativ zu ihrem Q, und 
homogen - kont inui er 1 ich uber 
alle Q relativ zu einem 
externen Bezugspunkt. 


Deformation wird als 
mat hemat ische Trans format ion 
("mapping") von Massepunkten 
aus den undef ormierten in die 


Deformation wird als 
e ne rge t i s che Zus t ands ande rung 
begriffen. Zwei verschiedene 
Punkte P und Q konnen nicht 



i 
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def ormierten Koordinaten 


auf einanderf alien , weil 


die 


begriffen. Urn zu verhindern, 


Arbeit sglei chung logarithmisch 


dass zwei Massepunkte bei der 


ist und dafur unendlich 


viel 


Transformation auf einander- 


Arbeit notwendig ware. 


Der 


fallen, ist der Apparat der 


Apparat 


der 


Kompatibilitatsgleichungen 


Kompatibilitatsgleichungen 


ist 


erforderlich. 


nicht erforderlich. 





Symbol terminologie 



T, M, F Tensoren f, m, r, s, v Vektoren 

I Identitatstensor n, t 

Einheitsvektoren 

T±j, F±j Tensorkomponenten r± 

Vektorkomponenten 

P, V, H, U Skalare, Zustands- S r - r- 

VektorgroSen 

Funktionen 

p, Q Punkte a, y, Q Winkel 



5 



10 



15 

Die Krafte, die die Oberflache aller P verschieben, setzen sich 
zusammen aus der Summe von externen Kraften f ex t, den vom System 
auf die Umgebung ausgeubten Kraften m^t, und den Oberflachen- 
Bindungskraf ten m*, die selbst keine Arbeit leisten (im Sinne 
20 einer Zwangskraft) , aber Dehnungskomponenten von f ex t/ sowie den 
Scherkomponenten m^yst) und f S (ext) es ermoglichen, Arbeit 
leisten zu konnen: 

fext + JHsyst + ™A = 0 • (Gl.01) 

25 

Neu ist an dieser Methode bei der Betrachtung von Spannung und 
Deformation 



/ 
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- die Einfuhrung eines fest definierten Bereiches (das 
System) , das von der Umgebung gedanklich geschieden ist 
und mit ihr in mechanische Wechselwirkung tritt, und 

5 dessen physikalische Eigenschaf ten maSstabunabhangig 

(unabhangig von der GroSe des Systems) sind, 

- die explizite Trennung von Kraften f ex t/ die von aussen 
kontrolliert von der Umgebung auf das System wirken, und 

10 Kraften m^st, die das System auf seine Umgebung ausubt, 

- die Einbeziehung der Existenz von Bindungskraf ten 
(permanent im Festkorper, nicht -permanent in einer 
Flussigkeit) , die die Freiheit der Reaktionen des Systems 

15 auf aussere Krafte erheblich einschranken und 

- die Erfassung der ^ni^-r^fc-i^Iien Wirkung von Scherkraf ten 
« • ,.***' .**■ ** 

(der von Scherkraft^ geleisteten Arbeit) und der 

Umwandlung der Scherkomponenten der Felder f ex t und nisyst in 
20 Normalkomponenten f fl (ext) und m^yst) parallel zur 

Hauptdehnungsrichtung . 

Die Erfindung sieht somit vor # ein Verfahren zur Ermittlung von 
Verformungen eines Korpers, wobei untersucht wird, wie die 

25 Einwirkung von Kraften den Korper verformt, so durchzuf iihren, 
dass fur einen Teil des Korpers oder seine Gesamtheit ermittelt 
wird, wie in Abhangigkeit von Materialeigenschaf ten, ausseren 
Grenzbedingungen und Vorhandensein von Bindungskraf ten zwischen 
System und Umgebung ein Kraf tegleichgewicht zwischen inneren 

30 und ausseren Kraften zu Stande kommt . Die Materialeigenschaf ten 
mussen als Funktion eines Ortes Q bekannt sein. Ebenso mussen 
die externen Grenzbedingungen (Eigenschaf ten des externen 
Kraftfeldes) bekannt sein. 
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Die gesamte von einem ausseren Vektorfeld geleistete Arbeit 
wird durch Aufteilung der von Normal- und Scherkomponente eines 
Vektors oder Vektorf eldes geleisteten Arbeit nach der Gleichung 



r * f = ^/|rxf| 2 +|r-f| 2 = |r||f | = const ( Gl . 1 8 ) 



ermittelt . 



Der Radius r ist ein mechanischer Hebel . Innerhalb eines 
10 Zeitraums, der verglichen mit der materialabhangigen 
Dif fusionsgeschwindigkeit bzw. Relaxationszeit kurz ist, ist 
das Verfahren auch auf Flussigkeiten und Gase anwendbar. 

Die Gesamtverschiebung wird nach einer logarithmischen 
15 Arbeitsglei chung ermittelt, die aus der dem Material 
entsprech^n^-- Zustandsgleichung herzuleiten ist ;^ im 
einf aclY;££SEr. -ist die Arbeitsgleichung 

jfdr = cjy = c\nr , (Gl.20) 

20 

die sich aus PV = nRT ableiten lasst. 



Nachfolgend wird die Erfindung anhand von zweidimensionalen 
25 Beispielf alien erlautert . Eine Erweiterung auf drei Dimensionen 
ist jedoch ohne Schwierigkeiten moglich. 

Weitere Vorteile, Besonderheiten und zweckmaSige 
Weiterbildungen der Erfindung ergeben sich aus den 
30 Unteranspruchen und der nachf olgenden Darstellung bevorzugter 
Ausfuhrungsbeispiele anhand der Zeichnungen: 



Von den Zeichnungen zeigen: 



WO 03/081484 PCT/DE03/01002 



Abb. 01 Oben: Bisherige Methode (Finite Element -Methode ) : ein 
Netz von Knotenpunkten wird vordefiniert und die Deformation am 
Punkt Q (im Zentrum) durch Anderung von Orient ierung und Lange 
5 aller Strecken bestimmt, deren eines Ende Q ist . Die jeweils 
ausseren Punkte sind dann Korrespondenzpunkte P. Die finite 
Element -Methode heiSt so, weil die Strecke zwischen zwei Knoten 
finit ist. Unten: Jeder beliebige Punkt kann als Ref erenzpunkt 
Q angesehen werden, alle Q zugeordneten Korrespondenzpunkte P 
10 bilden dann eine geschlossene Flache urn Q, und alle relevanten 
Variablen - Materialeigenschaf ten, Grenzbedingungen, Verformung 
- konnen als dif f erenzierbare Funktionen des Ortes f{Q) 
dargestellt werden. 

15 

Abb. 02 Kompressibilitatsdaten fur die Alkali -Elemente . Aussere 
Belastung in Vielfachen des inneren Drucks. Daten umgerechnet 
aus Bridgman (1958, p. 180) . 

20 Abb. 03 Beziehung von f op zu f d ev. Senkrechte Achse: 
VektorgroSe, horizontale Achse: Winkelabstand. Die linke und 
rechte Seite reprasentieren eine kontrahierende und eine 
extendierende Eigenrichtung (c und e) eines anisotropen Feldes. 
Nur in einem orthogonalen Feld sind die Eigenrichtungen 

2 5 orthogonal und gleichzeitig parallel zu den Hauptachsen eines 
Ellipsoids. Abb. 03 zeigt Beziehung der minimalen und maximalen 
ausseren Krafte relativ zum operativen Feld. Schwerkrafte sind 
nicht berucksichtigt . Umgebungsdruck und operatives Feld sind 
isotrop-hydrostatisch . 

30 

Abb. 04 (a) Deviatorisches Kraftfeld bei reiner Scherung. fi 
und f 2 sind senkrecht zueinander und parallel zu den 
Koordinaten und zu den Eigenrichtungen des Feldes. Dieselbe 
Abbildung reprasentiert das Verschiebungsf eld oder FlieSfeld, 
35 wenn die Pfeile fur Verschiebungen stehen. (b) Orthorhombisches 



WO 03/081484 




PCT/DE03/01002 



9 



gesamtes Kraftfeld, schematisch. m±: Materialvektor , f op : 
hydrostatischer operativer Druck, fn deviatorisches Kraftfeld. 

Abb. 05 KraftgrdSen entlang eines Einheitskorpers bei reiner 
5 Scherung. Durchgehende Linie: f n / gestrichelte Linie : f 8 . Ihre 
Summe ist 1 an alien Punkten. 



Abb. 06 Dilatanzef fekt einer Scherkraf t . Gebogene Linie: 
Oberflache des thermodynamischen Systems. Tangentialkraf tvektor 
10 (gestrichelt) wird Materialpunkt P von P nach P' verschieben; 
die Entfernung vom Ursprung vergrofiert sich. Punktkoordinate Xi 
wird von horizontaler Komponente von f auf xi' gelangt, x 2 wird 
von vertikaler Komponente von f auf x 2 ' verkurzt. X± : externe 
Koordinaten parallel zu den Eigenrichtungen. 

15 

Abb. 07 Aussere dicke Linie: Oberflache eines , B kxSs *5.k6 rp e r s\ . 
vor der Belastung. Innere dunne Linie: • Einv^.s'-v^Iu^ p&rr nachr 
hydrostatischer Belastung. aufgrund eines opea?^?5 V*?*s Feldes^? •••» 
Volumenschwund . Innere elliptische Linie: Formanderung des 

2 0 belasteten Korpers nur aufgrund der Normalkomponente des 

deviatorischen Feldes; Volumenkonstanz . Aussere elliptische 
Linie: finite Verformung nach Einbeziehung der von Scherkraf ten 
geleisteten Arbeit. Flache zwischen den elliptischen Linien: 
Volumenzuwachs aufgrund des Dilatanzef fekts . X±i aussere 
25 Koordinaten parallel zu den Eigenrichtungen. 

Abb. 08 Ausseres Kraftfeld fur einfache Scherung mit 
Einheitskorper. (a) Das Kraftfeld besteht nur aus Kraften 
parallel zu X x , (b) Der Einheitskorper ist im Raum fixiert und 

3 0 kann nicht rotieren. P: Ansatzpunkt eines Vektors, O: 

Koordinatenursprung, 0: Winkel zwischen Positionsvektor von P 
und Xi . 

Abb. 09 Krafte auf der Oberflache eines Volumenelements in 
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einem Kontinuum und einem Kraftfeld wie in Abb. 08 unterworfen. 
Durchschnittliche Kraft f av an ihrem Ansatzpunkt P K auf der 
Oberflache eines Einheitskdrpers bei einfacher Scherung. f av 
mit GroSe sin k zerlegt sich in Normalkomponente f n und 
Scherkomponente f 8 . Letztere ist ausgeglichen durch 
Oberflachenbindungskrafte. f 8(K > wird von alien Scherkraften an 
alien Punkten abgezogen, was sinistrale Scherkrafte zwischen P K 
und 0 zur Folge hat. X ist der Winkel einer 
Koordinatentransformation, die aufgrund des Ungleichgewichts 
der Scherkrafte zustande kommt . Die Gleichgewichtsbedingung 
J f x r d9 = 0 ist eingehalten. 

Abb. 10 Thermodynamisches System bzw. Volumenelement mit 
Einheitsmasse und elliptischer Form, die im Gleichgewicht mit 
einem monoklinen externen l^r^f-eXd ist. Spitzer Winkel 
zwischen Eigenrichtungen c ' -~ 

Abb. 11 VektorgroSen bei einfacher Scherung. (a) im 
untransformierten Zustand, (b) im transf ormierten Zustand. 
Durchgehende Linie: Normalkomponente fn(dev) (Gl.44, 
langgestrichelte Linie: | f d ev x r e ii I (Gl.42), Punktstrichlinie : 
f s (G1.45, siehe Text). J-R, -LP : Orientierungen senkrecht zu R- 
and P-Flache, vgl . Abb. 12. 

Abb. 12 Geometrie der einfachen Scherung. (a) Berechnete 
geometrische Eigenschaf ten: x it X±: Koordinaten, c: 
kontrahierende Eigenrichtung, e: extendierende Eigenrichtung, 
±R: Richtung vom Ursprung zum R-Punkt und senkrecht zur 
R(iedel) -Flache; JLP : Richtung vom Ursprung zum P-Punkt, 
senkrecht zur P-Flache. (b) Beobachtete geometrische 
Eigenschaften (Weber & Juckenack 1990) . Der einzige grofiere 
Unterschied zwischen (a) und (b) ist die Orientierung der 
Kornformfoliation. Es wird jedoch angenommen, dass die 
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extendierende Eigenrichtung mit der Gefuge -Orient ierung 
korreliert werden muss, die im allgemeinen etwas tiefer liegt. 

Abb. 13 (a) Verschiebungsfeld fur einfache Scherung. Obwohl das 
5 externe Kraftfeld keine Komponente parallel zu x 2 hat, ist eine 
solche Komponente im effektiven Feld f d ev und dem 
Verschiebungsfeld vorhanden. Eigenrichtungen und Koordinaten 
siehe oben. (b) Naturliches SC-Gefuge: Mylonit in unterer 
Grunschieferf azies von der Insubrischen Linie, Sesia-Zone, Val 

10 Strona, Italien. Lange Dimension des Photos ist parallel zur 
Gesamtfoliation, die als externes Bezugssystem dient . Die maSig 
nach rechts fallenden Scherflachen sind die C-Flachen. Der 
schwach nach links fallende stoffliche Lagenbau ist die S- 
Flache. Erhebliche sinistrale Gesamtkorper- Rotation zwischen C- 

15 Flachen-Diskontinuitaten, die selbst dextralen Schersinn haben. 
Gesamter Sche^ : inn ist dextral. (c) Viskose eiiifache Scherungf 
in subrez>v;. - ~.r ; hoidianf luss, Lipari/ Italien. Obere Lage .war 
weicher, b^3^^ r " -aiis schwarzem Glas mit Blasen und zeigt 
Schleppung. Untere Lage war steif, besteht aus teilweise 

20 kristallisiertem Material und reagierte bruchhaft. Schleppung 
in oberer Lage und Kluf torientierung in unterer Lage zeigen 
dextrale Scherung an. 

Abb. 14. Oben: rechteckiger Festkorper zwischen zwei 

25 Stempeln mit den Orten, fur die die Grenzbedingungen unten 
dargestellt sind. Mitte: Grenzbedingungen als Funktion des 
Ortes im Korper. Kreise stehen fur (von innen nach aussen) : 
operatives Feld (isotrop) ; Normal komponente des deviatorischen 
Feldes; vom System beigesteuerter Scher-Dehnungsef f ekt ; 
3 0 Volumenef fekt aufgrund der externen Grenzbedingungen. Letzteres 
ist in (1) wegen des Reibungsef f ekts am Werkstuck-Stempel- 
Kontakt einwarts gerichtet, in (2) wegen des von der Umgebung 
beigesteuerten Scher-Dehnungsef f ekts auswarts gerichtet, und 
Null in (3) fur Punkte auf der Grenzflache zwischen Festkorper 
35 und freiem Raum. Unten: Darstellung des erwarteten Effekts; 
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weisse Kreise: unbelasteter Zustand, schwarze Ellipsen: 
belasteter Zustand. 



Abb. 15. Dimensionen 10x2 Entf ernungseinheiten. Wellige 

5 Struktur entlang der horizontalen Rander im unteren Bild sind 
mathematische Artefakte aufgrund Summierung uber nicht 
unendliche Koef f izienten n. 



Abb. 16. Dimensionen 10x5 Entf ernungseinheiten. Wellige 

10 Struktur entlang der horizontalen Rander im unteren Bild sind 
mathematische Artefakte aufgrund Summierung uber nicht 
unendliche Koef f izienten n. 

Abb. 17 Dimensionen 10x10 Entf ernungseinheiten. 

15 

Abb. 18 Dimensionen 10x30 Entf ernungseinheiten. 



Modellrechnungen sind hier aus Grunden der Vereinf achung auf 

zweidimensionale Beispiele beschrankt. In der vorliegenden 

Anmeldung sind die Gleichungen in der Regel in Vektornotation 
2 5 und in algebraischer Form gegeben, z.B. 

J f x r d9 = J sinG cosG d9 



Die Erfindung umfasst insbesondere folgende vier Teile: 
3 0 - der Definition eines maSstabunabhangigen Einheitsvolumens 

mit einer den Gleichgewichtsbedingungen entsprechenden 
Form, 

- der Definition zweier Kraftfelder, die miteinander im 
Gleichgewicht stehen, 



[ 
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- der Aufstellung einer materialunabhangigen 
Zustandsglei chung fur Festkorper, 

- der Umformung der Zustandsgleichung von Skalar- in 
Vektorf orm. 

5 

Definition des Systems 

Zweck des Systems ist die Bereitstellung eines Ref erenzkorpers 
10 im unbelasteten Zustand, der die Material- und andere 
mechanische Eigenschaf ten reprasentiert . Auf dieses System wird 
bei Belastung Arbeit geleistet, die maSstabunabhangig ist in 
dem Sinn, dass zum Beispiel auf ein kleines und ein groSes 
System bei gleicher Verformung relativ zu seiner Masse gleich 
15 viel Arbeit geleistet wird. Damit wird gewahrleistet , dass man 
zwar mit einem finiten Volumen oder einem Radius reefed kann r ,;i 
die tatsachliche GroSe des Systems jedoch nicht re l&v.? r.r. - i&tv <35*S 

Jegliche Masse ist mit einem Potential assoziiert. Im hier 
20 betrachteten Fall wird die innere Energie eines 
thermodynamischen Systems als sein Potential betrachtet. Sie 
ist in einem gegebenen energetischen Zustand proportional zur 
Masse, dU oc dn, wobei n die Zahl der Mole ist (die 
thermodynamische Masse ist dimensionslos) . Die Energiedichte 
25 pro Masse ist dann dU/dV = P, in integrierter Form U/V = P, 
wobei die Masse proportional sowohl zu U als auch zu V ist. U/V 
= P ist die Energiedichte, oder das thermodynamische Potential 
pro Einheitsmasse . Letztere ist dann finit. 

3 0 Im Standardzustand C7 0 gilt PV = nRT. Eine inf initesimale 
Zustandsanderung erfordert die zusatzliche Energie dl7. Dieses 
dU ist bei konstantem n von Natur her etwas ganz anderes als 
dU oc dn oben. Somit ist ein System von finiter Masse und 
finitem Volumen, assoziiert mit seinem Eigenpotential , eine 
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Voraussetzung, um inf initesimale Zustandsanderungen zu 
betrachten. 

Um auf einen Punkt P 0 wirkende Arbeit definieren zu konnen, 
5 muss ein Ref erenzpunkt Q * P 0 gewahlt werden. Die Wahl ist 
beliebig und kann sich aus der Natur des Problems ergeben 
(Kellogg 1929, p. 53, 63; Boas 1983), Die Basisentf ernung Q-»P 0 
ist die Nullpotentialentfernung r und wird hier interpretiert 
als der Radius des thermodynamischen Systems, im anisotropen 

10 Fall: als die Entfernung eines jeden Punktes auf der Oberflache 
des thermodynamischen Systems zu dessen Massezentrum. Da die 
Masse gleichmaSig in der vom thermodynamischen System 
eingenommenen Region verteilt ist, stellt sie eine verteilte 
Quelle im Sinn der Potentialtheorie dar (Kellogg 1929, p. 156). 

15 Die Potentiale von verteilten Quellen sind allgemein 
logarithmisch. Die Nullpotential-^n^f &^nung kann dann jeden- 
f initen Wert (ausser Null) habeas*. ^ur;i. aber per Konvention <S£t~ 
auf 1 gesetzt . 

2 0 In der Thermodynamik werden die Eigenschaf ten der Materialien 

in einem Standardzustand (P*, V* , T* ) als gegeben, der dann 
auch als Nullpotentialzustand angenommen werden kann. Eine 
Druckanderung auf absolut 0 oder oo erfordert dann unendlich 
viel Arbeit. Der Entf ernungsterm r wird hier interpretiert als 
25 der Radius r 0 des thermodynamischen Systems im Standardzustand. 
Er ist das eindimensionale Aquivalent des Volumens der Masse n 
in der Zustandsgleichung PV - nRT. Fur Festkorper, die im 
Vakuum ein finites Molvolumen und eine Dichte > 0 haben, ist 
der Nullpotentialzustand dann definiert als P* = Pint- Eine 

3 0 Veranderung der GroSe von r 0 erfordert dann eine 

Zustandsanderung . In erster Naherung wird r 0 zunachst als 
invariant mit Bezug zur Richtung angesehen, d.h. es wird 
angenommen, dass das thermodynamische System Kugelform und die 
Oberflache pro Volumen einen Minimalwert hat. Die 
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Grenzbedingungen wie z . B . 
fuhren, dass eine andere Form 



5 Im isotropen Fall brauchen nur Normalkrafte betrachtet werden, 
so dass der Divergenzsatz 

Jf ndA = JV f dV = K (G1.02) 
die Situation vollstandig beschreibt. n ist der Einheitsvektor 
senkrecht zu A. Kraft wird als Energiefluss f = & i dU/dx i 

10 interpretiert . A ist eine geschlossene Flache, die das 
Volumenelement JdV enthalt; die Krafte werden von der Masse im 
System auf die Umgebung ausgeubt, als auch umgekehrt, d.h. es 
sind zwei Kraft f elder, die miteinander in Wechselwirkung treten 
und jedes fur sich Gl.02 erfullen (wobei k jeweils 

15 entgegengesetztes Vorzeichen hat) . Unter gleichen aussereri 
Bedingungen mu>V\s o^r.Ai kleine und eine groSe Masse relativ zur 
betrachteten gleiche Arbeit auf ihre Umgebung leisten, 

die Betrachtung muss also mafistabunabhangig sein. Das 
thermodynamische System stellt daher eine Quelle von Kraften 

2 0 dar. Die Divergenz des vom System ausgeubten Kraft feldes (und 
daher auch die Divergenz des von aussen auf das System 
ausgeubten Feldes) ist proportional zur Masse, also ist die RHS 
eine lineare Funktion von Vj daher ist div f eine Konstante . 
Daraus f olgt : bei gegebenen externen Bedingungen ist div f 

2 5 maSstabunabhangig. Wenn der betrachtete MaSstab, also V, 
verandert wird, kann, da die Relation von V oc r 3 zu A oc r 2 
nicht linear ist, das Verhaltnis von A zu |f| nicht konstant 
sein; |f| muss notwendigerweise eine lineare Funktion des 
MaSstabs sein, |f| oc r. div f ist die Spur des unten 

30 definierten Tensors F. Sowohl aus Gl . 02 und G1.06 folgt 

M = divf = const (G1.03) 
M 

unter konstanten externen Bedingungen, wobei r der Radius des 
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Systems und r ein Mass fur den betrachteten MaSstab ist. 

Diese Schlussfolgerungen werden nicht verandert durch weitere 
Entwicklungen des Divergenzkonzeptes , die erst weiter unten 
5 erklart werden konnen, wenn der Kontext deutlicher ist- Gl.03 
ist vereinbar mit den Eigenschaf ten eines thermodynamischen 
Kontinuums, mit den raumlichen Eigenschaf ten eines 
thermodynamischen Systems, das mit seiner Umgebung in 
Wechselwirkung tritt, und mit der Tatsache, dass das System ein 
10 verteiltes Potential darstellt (Kellogg 1929, p. 156). 

Ein thermodynamisches System wird durch eine gegebene Menge an 
Masse definiert. Sein Ort im Raum ist durch sein Massezentrum 
(oder seinen Schwerpunkt) in einem externen Koordinatensystem 

15 X ± gegeben. Die Oberf lachenpunkte des Systems werden mit Hilfe 
eines eingeschriebenen, inneren Koordinatensystems Xi 
A^viVttelt, dessen Ursprung mit dem Massezentrum des gewahlten 
Systems zusammenf allt . Die physikalischen Bedingungen an den 
Punkten P in xi sind Funktionen von Q in Xi im Sinne eines 

2 0 Vektorfeldes, T(Q)v = x, wobei T ein Tensor und eine 
Ortsfunktion von Q ist; v ist der Ortsvektor eines in der Nahe 
liegenden Punktes P relativ zu Q, und x ist ein Vektor am Ort P 
als Funktion von T und v. 

2 5 Die Beziehung zwischen einem Oberf lachenpunkt P und einer 
bestimmten Richtung in xi ist eindeutig. Ein anderer Vektor 
durch P kann nur ein Teil eines anderen Systems mit seinem 
eigenen Ursprung und seinen eigenen Koordinaten sein, und es 
mogen andere Grenzbedingungen gelten, falls die Gradienten in 

30 Xi nicht null sind. Er ist daher zunachst nicht von Interesse. 

Ableitungr des ausseren und des Materialkraftfeldes 



3 5 Eine Definition eines Tensors ist die Ableitung eines 



10 
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Vektorfeldes nach den Koordinaten, oder die zweite Ableitung 
eines Skalarfeldes nach den Koordinaten (Malvern 1969, p. 57). 
Energie oder Arbeit sind solche Skalare. Es wird angenommen, 
dass ein externes Potential U ex t exist iert, und die Material - 
Enthalpie sei E, dann ist 

™~L=f. — = (G1.04) 

8x l J ' 

^. = ^L = efF . e , _^ = ^L = e( Me, (G1.05) 



_ _ _ F • ej 

dxfixj dxj dxgdxj dx ; 



(G1.06) 



wobei r der Radius bzw. der Ortsvektor eines Punktes auf der 
Oberflache des Systems, f das aussere Kraftfeld, m das 
Materialkraftfeld des Systems (einfach: Materialf eld) sind, und 
P und M sind die Eigenschaf tstensoren der beiden Vektorf elder f 
und m. F represent iert die ausseren Grenzbedingungr;£ = i d m 
stellt die Materialeigenschaf ten dar, die daher auci ■ 
15 Art innerer Grenzbedingung verstanden we r den konnen. Aur£; das 
System wirkende Krafte werden als Energiefluss verstanden, der 
den inneren Energiezustand verandert, nicht aber als externe 
Beschleunigung eines Punktes. 



20 Da ausserlich unausgeglichene Krafte keine Deformation bewirken 
konnen, werden sie ignoriert . Ausserdem werden Massekrafte 
ignoriert, die mit der tragen Masse in Wechselwirkung treten, 
aber nicht mit der thermodynamischen Masse- Die innere 
Gleichgewichtsbedingung ist dann gegeben durch 

25 f + m = 0 (G1.07) 

fur jeden Punkt P auf der Oberflache des Systems, oder, wenn f 
und m als Funktionen der Richtung G verstanden werden, 

jfde = -jmdQ (G1.08) 
als Summe uber die Oberflache des Systems in 2D. 



30 
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Das aussere rotationale Gleichgewicht ist per definitionem 
eingehalten, da System und Umgebung fest miteinander verbunden 
sind; der Ungleichgewichtsf all kann in einem elastischen Medium 
nicht eintreten, so lange keine Bindungen gebrochen werden. Die 
Bedingung 

Jf xrd9 = 0 (G1.09) 

ist daher keine Gleichgewichtsbedingung, sondern lediglich die 
Feststellung der Orthogonal it at fur den Fall, dass das 
rotationale Gleichgewicht eingehalten wird ohne zusatzliche 
Hilfe durch Bindungskraf te n^. Die vollstandige Bedingung fur 
das rotationale Gleichgewicht ist 

J(f xr - m s xr)t/9 = 0 (G1.10) 

wobei nig 1 r, ihre GroSe ist invariant, 

ffcJ = Q (Ol.ll) 

ae 

und ms x r hat dasselbe Vorzeiche^ * .. r+bX*& Oberf l^chenpurikfteri 
P. Bindungskraf te sind weder ato^rliche Krafte hbch 
Materialkrafte; sie exist ieren unabhangig von den 
Materialeigenschaften. Sie treten nicht in Erscheinung in 
Wechselwirkung zu einwarts-gerichteten (komprimierenden) 
Normalkomponenten von f, sondern nur von auswartsgerichteten 
(dilatierenden) Komponenten sowie allgemein mit Scherkraf ten; 
ihre Existenz folgt aus der Vorbedingung, dass Gleichgewicht 
existieren muss, und sie - falls notig - das Drehmoment von f 
ausgleichen. Das Thema wird spater genauer illustriert. 

G1.07 und G1.08 sind wesentlich fur das Verfahren, weil sie 
aussere und Materialkrafte klar unterscheiden, d.h. 
Materialkrafte als eigenstandige Krafte erkennen. Sie werden 
vom System auf die Umgebung ausgeubt aufgrund einer 
Zustandsanderung im System aufgrund der Einwirkung ausserer 
Krafte. Das System stellt daher ein Potential dar. 
Gleichgewicht zwischen System und Umgebung bedeutet unter 
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isotropen Bedingungen, dass 

div f + div m = 0 ; (Gl . 12 ) 

G1.12 besteht daher aus zwei Poisson-Gleichungen. 

5 Spannung im Sinne dieses Antrags ist ein Kraf tvektorf eld f total/ 
das aus der Wechselwirkung zweier unabhangiger Kraftfelder f 
und in, dem ausseren und dem Materialkraf tf eld, ihren jeweiligen 
Grenzbedingungen im Gleichgewichtszustand, plus der 
Bindungskrafte falls erforderlich besteht (vgl.Gl.01), 
10 ftota, =/(f, m, m s ). (G1.13) 



Aufstellung einer materialunabhangigen Zustandsgleichung fur 
Festkorper 

15 

Das ideale Gasge^e^r^ i?V = nRT, ignoriert die atomistische 
Struktur von und betrachtet nur makroskopisches 

Gesamtverhalten. DW Sbylesche Gesetz kann als Materialgesetz 
fur ein ideales Gas betrachtet werden, da es ein bestimmtes 

20 Verhalten des Gases aufgrund einer Zustandsanderung vorhersagt . 
Es bringt die innere Energie eines Systems, seine Masse und 
sein Volumen miteinander in Beziehung und beschreibt daher den 
ihm eigenen energetischen Zustand vollstandig- Im Prinzip muss 
eine solche Zustandsgleichung auch fur Festkorper existieren. 

25 Es muss daher moglich sein, einen idealen Festkorper zu 
def inieren. 



Als allgemeinste Form der Zustandsgleichung fur einen 
Festkorper wird 

30 P k V = z (01.14) 

gewahlt, wobei 

lnK!° Hd 



k= ™' • (G1.15) 



Die Formel besagt, dass alle Festkorper dieselbe 
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Kompressibilitat dV/dP aufweisen, wenn das molare Volumen 
berucksichtigt wird. z = f(P) ist eine Zahl, die fur einen 
gegebenen Zustand charakteristisch ist; die Funktion selbst ist 
nicht bekannt, wird aber mit Hilfe der Birch- Murnaghan- 
Gleichung erfolgreich modelliert. Diese wird bei der 
Model lierung von Materialverhalten unter hohen Drucken 
angewendet, z.B. der Vorhersage der elastischen Eigenschaf ten 
des Erdkerns (Jeanloz & Knittle 1986) . 



Es ist bekannt, Alkalimetalle auf 100 kb zu komprimieren. Sie 
folgen in erster Annaherung alle derselben Funktion, wenn die 
ausserlich angelegte Belastung in Vielfachen des inneren Drucks 
des jeweiligen Materials ausgedruckt wird (Abb. 02) . Es ist 
daher gerechtf ertigt , den inneren Druck allgemein als Standard 
fur Materialverhalten zu benutzen. Gl . 14 ist isotrop, daher ist 
k -ichtungsunabhangig, und ausserdem nur wichtig fur das 
f/'.v„/;.;eren von echten Festkorpern. Es genugt, wenn k = 1 
■'■gm*®«<zt wird; damit ist z = const. Somit sind die 
thermodynamischen Eigenschaf ten des idealen Festkorpers 
identisch mit denen eines idealen Gases, und das Boylesche 
Gesetz ist anwendbar. 



Umformung der Zustandsgleichung von Skalar- in Vektorform 

Das fundamentale Materialgesetz in der Thermodynamik ist das 
Boyle -Mariottesche Gesetz, 

PV = const (G1.16) 
Beide GroSen links konnen neu interpretiert werden, ohne dass 
die Einheiten selbst sich andern (es bleibt das Joule) ; das 
Ergebnis ist eine Vektor-Form des Boyleschen Gesetzes, 

r*f = const (G1.17) 
wobei r der Ortsvektor eines Punktes P auf der Oberflache des 
Systems relativ zu seinem Massezentrum ist. Das Sternprodukt 
wird interpretiert als 
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r*f = 7|rxf| 2 +|r.f| 2 = |r||f| = const ; (Gl . 18 ) 

r*f beschreibt die Arbeit, die von der Normal- als auch der 
Scherkomponente einer externen Kraft f am Punkt P geleistet 
wird, der den Ortsvektor r hat und Teil der Oberflache eines 
5 thermodynamischen Systems ist. G1.18 ist eine bekannte 
Identitat. Wenn r ein Einheitsvektor ist, gilt r*f =|f|; aber 
das Verhaltnis von |r| zu |f| hat nun die Freiheit, als 
Funktion der Richtung variabel zu sein, auch wenn G1.18 
eingehalten wird, wenn namlich |f| und |r| sich zueinander 
10 umgekehrt reziprok verhalten. 

Das Kraft f eld f total/ das sich beim Ubergang vom unbelasteten in 
den belasteten Zustand entwickelt, kann in eine isotrope 
Komponente, das operative Kraftfeld f oP / und die deviatprische 
15 Komponente f de v zerlegt werden, wobei fop die Zustandsande^m^c; 
reprasentiert . Sowohl f-r als auch |f x r|integriert ub^r <£2£& 
Oberflache des Systems summieren sich auf Null. Durch f op • 
ein Referenzzustand definiert, der dazu dient, den Richtungen 
des deviatorischen Feldes Vorzeichen zuzuordnen (Abb. 03) 

20 

Im Gegensatz zu Gasen und Flussigkeiten konnen Festkorper 
Scherkrafte aushalten, und der energetische Zustand eines 
Systems wird von Scherkraften und Normalkraf ten gleichermaSen 
beeinflusst. Eine Normalkraf t wird je nach Vorzeichen eine 

25 Dilation oder Kontraktion des Systems verursachen. Das 
Vorzeichen einer Scherkraft hingegen zeigt nur ihren Drehsinn 
an, aber es hat keine Aussage uber einen Volumenef f ekt . Die 
Wirkung der Scherkrafte auf das Volumen kann jedoch nur 
dilatierend sein. Beispiel : ein deviatorisches Feld f dev kann 

3 0 weiter aufgeteilt werden in ein nur aus Norma lkomponen ten f n 
und ein nur aus Scherkomponenten f s bestehendes Feld. Ein 
kugelformiges Volumenelement wird den Scherkraften f a eines 
Kraftfeldes mit orthogonalen Eigenrichtungen unterworfen; der 
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Ursprung des Koordinatensystems sei am Massezentrum. Auf der 
Systemoberf lache setzen die f s auf den Punkten P mit 
Positionsvektor r an und verschieben sie parallel zur Richtung 
von f s . Sie werden die Positionen P' mit Positionsvektor r 1 
5 haben. Da r und v = P -> P' senkrecht zueinander stehen, ist 
immer |r'| > |r|. Daher hat die von einer Scherkraft geleistete 
Arbeit einen Volumenef f ekt , der immer dilatierend ist. Die 
dilationale Komponente summiert sich von Null am Ansatzpunkt 
der kontrahierenden Eigenrichtung in Richtung auf den 

10 Ansatzpunkt der extendierenden Eigenrichtung, wo sie ihr 
Maximum erreicht. Diese Wirkung gilt sowohl fur Scherkrafte, 
die von der Umgebung auf das System ausgeubt werden, als auch 
umgekehrt fur die Krafte, die vom System auf die Umgebung 
ausgeubt werden. Der gemeinsame Effekt der Scherkrafte 

15 resultiert in zwei zusatzlichen normal -dehnenden Kraftfeld- 
Komponenten f 8( ext) und m s(sy8 t) -parallax '%Hr' extendierenden 
Eigenrichtung. Im Endergebnis wird a^;^*'.'. %r . Sxistlnz. der 
Scherkrafte eine dilatierende (also normals) Wirkung parallel 
zur extendierenden Eigenrichtung erzielt, so dass 

2 0 f total = f op + fdev + fs(ext) + ™s(syst) (Gl . 1 9 ) 

wobei fext = fop + fdev/ F op = cl , und det F dev = ± 1- Ein Beispiel 
wird unten gegeben. Die Tensoren sind Feldeigenschaf tstensoren 
wie in der Definition eines Vektorfeldes Tv = x in der linearen 
Algebra . 

25 

Wenn das gesamte Kraftfeld f total bekannt ist, kann die Arbeit 
analog zur PdV-Arbeit berechnet werden, 

J/^ = cJy = clnr (G1.20) 

Diese Gleichung, uber die Oberf lache eines thermodynamischen 

3 0 Systems integriert, liefert dasselbe Ergebnis wie die einfache 

Berechnung der PdV-Arbeit, wenn die Belastung isotrop und 
komprimierend ist. 
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Die Beziehung einer Materialdistanz, also r 0 , zu einer 
kolinearen Kraft ist durch Differentiation von G1.18 gegeben, 

r o df + f o dr = 0 (G1.21) 
Wenn Gl.21 durch r 0 dividiert, re-integriert und durch f 0 
dividiert wird f folgt 




(G1.22) 



in vollstandiger Analogie zu In (Vi/V 0 ) = - AP/P 0 in der 
isotropen Thermodynamic Gl.22 stellt die Gleichung dar, die 
Ursache und Wirkung miteinander verbindet; rait ihrer Hilfe wird 
das Verschiebungsf eld Ar = s aus dem Kraftfeld ftotai 
ausgerechnet, die somit identische geometrische Eigenschaf ten 
haben . 



Ein Festkorper hat im Vakuum sein ideales Volumen V 0 mit 
Einheitsradius r 0 . D^i : v ^>^:v > Druck P int eihes Festkorpers liegt 
in der GroSe von li-^JOl SK?a. Analog zum inneren Druck in 
Vektorform ist die innere oder Materialkraf t mo # die eine Kraft 
mit Grofie >0, aber im unbelasteten Zustand innerlich 
ausgeglichen ist. Zu Model lierungszwecken ist sie im Folgenden 
auf Einheitsgrofie gesetzt. Die Grofie der ausseren effektiven 
Kraft muss mit Bezug zur Grofie der Materialkraf t mo normiert 
werden, woraus Gl.07 vollstandig lautet : 



-Am = M CTt =i^ ] ^ (G1.23) 

l m o| l f ext| 

wobei der erste Term rechts die Grofie von f ex t in Mehrfachen des 
inneren Drucks des Materials, der zweite ein Einheitsvektor mit 
der Richtung von f ex t ist. Diese Normierung ist in Gedanken 
immer angewandt; im Folgenden werden der Einfachheit halber 
Einheitsgrofien benutzt, und Af wird als f bezeichnet. 



Die Gleichgewichtsbedingung fur das rotationale Moment ist 
abhangig von den Eigenschaf ten des Kraft feldes und der Form des 
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Systems. Das Gleichgewicht ist immer erhalten, weil System und 
Umgebung fest miteinander verbunden sind. Die Form des Systems 
wird bestimmt von den Eigenschaf ten des Kraftfeldes, den 
Mater ialeigenschaf ten und der Bedingung, dass das A/V- 
5 Verhaltnis bei gegebener Masse minimiert sein soil. 

Wenn die Wirkung einer progressiven Verformung in 2D gezeigt 
werden soli, wird zumeist ein Kreis durch ein Verschiebungsf eld 
in eine Ellipse verwandelt. Der Einheitskreis im unverformten 

10 Zustand ist sicher ein adaquates Konzept fur die Demonstration; 
aber die Form des thermodynamischen Systems (oder 
Volumenelements) muss die Gleichgewichtsbedingungen einhalten. 
Insbesondere kann das System im undef ormierten Zustand eine 
elliptische Form unter einem bestimmten Satz von Grenz- 

15 bedingungen haben, in einem anderen Zustand kann es ein Kreis 
sein. Pie beiden Konzepte - eine kreisformige Vorzeichnung im 

un-<J&f r-r / ^:ten Zustand gg. die mechanisch aktive Form des^ 

* ■ -J*- ■ 

S#sk«&?^ durfen nicht verwechselt werden. 

2 0 Eine Funktion des Verschiebungsf eldes ist der Strain-Tensor. 
Wenn die Hauptachsen des Strain- Ellipsoids wahrend einer 
progressiven Deformation/ nicht rotieren, mussen sie mit den 
charakteristischen Richtungen des Verschiebungsf elds (und daher 
auch des Kraftfeldes) zusammenf alien. Es kann daraus 

2 5 geschlossen werden, dass die Eigenrichtungen des Kraftfeldes 

und daher auch das Verschiebungsf eld orthogonal sind. Aus 
mathematischer Sicht ist dies jedoch ein sehr spezieller Fall, 
und keineswegs eine Vorbedingung . Der Strain-Tensor ist daher 
nicht hilfreich beim Verstandnis von Deformation. Es ist daher 

3 0 notwendig, (a) die Form des Systems aus den 

Gleichgewichtsbedingungen zu schliessen, und (b) die 
Eigenrichtungen des Kraftfeldes zu finden. Sie sind definiert 
als die Richtungen, entlang derer das Spannungsf eld f keine 
Scherkomponenten hat . 



35 
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Form des Volumenelements 

G1.18 liefert Grenzbedingungen fur die Form des Systems. Die 
5 einfachste Form ist die Kugel - sie ist isotrop, minimiert das 
A/V-Verhaltnis, und fur jede Punktquelle am Massezentrum ist 
ihre Oberflache eine Aquipotentialf lache . Fur isotrope externe 
Verhaltnisse und ein isotropes Material ist die Kugelform fur 
das System die naturlichste Wahl, weil keine Scherkrafte 
10 existieren. Fur andere Verf ormungstypen schreibt G1.18 vor ; 
dass Radius und Normalkraft entlang der kontrahierenden 
Eigenrichtung c und entlang der extendi erenden Eigenrichtung e 
identische GroSe haben. Die Wirbelung des deviatorischen Feldes 
fur reine Scherung (Abb. 04) 



15 




ist Null; daher fallen die Eigenrichtungen des resultierenden 
Feldes mit denen von F dev zusammen. Die Bedingungen tr T = 0 und 
det T = ± 1 sind beide Konservierungsbedingungen und mussen 
gemeinsam eingehalten sein, also |Fnl = I F 22 1 = 1. Da F de v 
20 orthogonal ist, gibt es keine Beschrankungen aus den 
Gleichgewichtsbedingungen heraus fur die Form des Systems, das 
daher immer noch kugelf drmig ist. 

Die GroSen der Normal- und der Scherkomponenten der 
25 deviatorischen Krafte sind in Abb. 05 gezeigt. Unter den 
gegebenen Grenzbedingungen haben alle deviatorischen Krafte auf 
der Oberflache die Komponenten f = [cos9, -sin9] , 9 ist dabei 
der Winkel zwischen der Orientierung von r und der x x - 
Koordinate. Die externe Gleichgewichtsbedingung fur das 
30 Drehmoment fur den gesamten Korper ist dann 

jf'ff xr J9 = 2j^ n cose sin0rf9 = O, (Gl.25) 
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und die Oberf lachenbindungskraf te ma sind nicht aktiviert. 
Kinematik der reinen Scherung 

5 

mo ist die Materialkraf t , die mit einem Vakuum im Gleichgewicht 
ist. Fur einen Korper mit Einheitsradius in einem orthogonalen 
deviatorischen Feld sind die Komponenten des Posit ionsvektors r 
eines Oberf lachenpunktes P durch [r x r 2 ] = [cos9 sinB] gegeben; 
10 der operative Kraftvektor ist eine Normalkraft mit den 
Komponenten f op proportional zu [-cos0 -sin0] ; die 
deviatorische Kraft f dev ist proportional zu [cosG -sin0] ; und 
der radius- normal e Tangens-Einheitsvektor t = [sin0 -cosG] . 

15 Auf P wirkende externe Krafte sind die isotrope operative Kraft. 
f op und die deviatorische Kpmponente f dev < ^; : ;i?;eht nur, aus 

Normalkraf ten. f dev kann we iter . aufgete; l^i : -nv^den ^ in < eihe 
Normalkomponente f n (dev) und eine Tangentialkomponente f S ( deV ) • Die 
Summe aller Normal komponenten ist 

20 2f n = + (r-f^Jr = 1 + (cos 2 9 - sin 2 0) (Gl . 26) 

Der entsprechende einwarts gerichtete Verschiebungsvektor s n 
parallel zu r kann durch Gl.22 gefunden werden. 

Das Kreuzprodukt 

25 | r x | = If^l = 2 cos0 sin0 (Gl . 27 ) 

ist gleichwertig mit f • t. Eine Scherkraft bewirkt, dass P von 
der kontrahierenden Eigenrichtung weg- und zur extendierenden 
Eigenrichtung hingeschoben wird (Abb. 06) . Die Komponenten von 
fs(dev) bei P rufen im System Materialkraf te gleicher Grofie und 
30 Richtung, aber entgegengesetzter Ausrichtung hervor, so dass 
das Drehmoment bei P ausgeglichen ist. Die Wirkung von f 8 ist, 
dass P nach P» verschoben wird (Abb. 06); da jedoch f s senkrecht 
zu r steht, ist die Entfernung r 1 vom Ursprung nach P' stets 
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langer als die Entfernung vom Ursprung nach P (Abb. 06). Also 
wird eine Scherkraft stets eine lokale Langung proportional zu 
ihrer eigenen Grofie verursachen. Dieser Effekt ist additiv von 
Punkt zu Punkt. Die beobachtete Dilation bei P' ist daher die 
5 Summe aller Dilationen von der kontrahierenden Eigenrichtung 
bis P, 

J( f dcv • *) r dQ = 2 £ sin9cos9>/cos 2 6 + sin 2 9 dQ (G1.28) 
wobei a die Winkelentf ernung von der kontrahierenden 
Eigenrichtung c bei x 2 bis P ist. Fur a = */ 2 ist Jf-t dQ = 1 . An 
10 der unteren Integrationsgrenze gibt es naturlicherweise keine 
dilatorische Wirkung; das Integral erreicht seinen hochsten 
Wert an der extendierenden Eigenrichtung e, obwohl an diesem 
Punkt keine Scherkrafte vorhanden sind. 

In einem allseits verbunder."-; Kontinuum zahlt das Ergebnis von 
G1.28 doppelt. Unter dev xS<-&, ein Korper sei hydrostatiscfh 
belastet, so dass gg****- operative Kraftfeld (das die 
Zustandsanderung bewirkt) die Grofie -1, und das deviatorische 
Feld die Grofie Null hat. Entlang von x x und x 2 wird der Radius 
verkurzt. Dann werden die Grenzbedingungen entlang von x x 
gelockert; es wird angenonimen, dass die Umgebung immobil ist. 
Im Versuch, die gespeicherte Energie zu reduzieren, expandiert 
das Material von selbst parallel zu Xi, aber nichts geschieht 
parallel zu x 2 . Die GroSe des operativen Feldes reduziert sich 
auf -X, und ein deviatorisches Feld mit der GroSe ± entsteht . 
Da die Ausdehnung nicht isotrop ist, ubt der Korper Scherkrafte 
auf seine Umgebung aus; rechts von der kontrahierenden 
Eigenrichtung wird ein linksdrehendes Kraftepaar an der 
Zwischenf lache von System und Umgebung oberhalb der 
extendierenden Eigenrichtung (im 1 . Quadranten) beobachtet, 
darunter aber ein rechtsdrehendes (im 2 .Quadranten) . Also wirkt 
das System auf seine Umgebung, wobei der energetische Zustand 
des Systems von -1 auf reduziert wird, bis das System seinen 
energetischen Minimal zustand erreicht hat. Allerdings ist zu 
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bemerken, dass das Vorzeichen der Scherkraf tepaare 
entgegengesetzt zu dem ist, was man bei einer reinen Scherung 
im Kontinuum normal erweise erwarten wurde. Wenn jedoch die 
Umgebung zusammen mit dem System verformt wird, wird sie erst 
5 mit dem System in x ± expandieren und dann zusatzliche 
Scherkraf te auf das System ausiiben, was im l.Quadranten ein 
rechtsdrehendes Kraftepaar und im 2.Quadranten ein 
linksdrehendes Kraftepaar bewirkt, was wiederum zu weiterer 
Ausdehnung in Xi f uhrt . Diese Scherarbeit ist von der Umgebung 

10 geleistet und fuhrt zu einem Anstieg des energetischen Zustands 
von -M auf -1. Also ist ein Teil der Ausdehnung vom Schereffekt 
des Systems, der zweite Teil vom Schereffekt der Umgebung 
beigesteuert , und beide sind proportional zur GroSe des 
deviatorischen Feldes. Im Endergebnis setzt also eine 

15 dilatorische Kraft der doppelten GroSe der integrierten 
ausseren Scherkraf te auf den Oberf lachenpunkt unter ;.. der 
extendx-it'^ii^^ Eigenrichtung an, der daher vom Massezentrum des 
SysteuiK •;r^rc|e.schoben wird. 

20 

Volumeneffekt der deviatorischen Belastung 

Es wird ein numerisches Beispiel ausgehend von G1.22 gegeben. 
r Q = 1, f 0 = 1/ und Volumen V 0 /n = 1. Nach Aufladung zur 
25 operativen Feldstarke -1/5 sind die GrdSen entlang x 2 : Af = fop 
+ -fn(dev) = -2/5, r 2 = 0,670, entlang x x : f op + ^n(dev) = -1/5 + 1/5 
= 0, r x = 1,000. Jedoch f 8( int) + ^s(ext) = 2/5, daher r x = 1,492. 
Das Volumen der resultierenden Ellipse ist V/n = r x r 2 = 1, d.h. 
das Volumen ist dasselbe wie im unbelasteten Zustand (Abb. 07) . 

30 

Model lrechnung fur einfache Scherung 



Die Eigenrichtungen 



WO 03/081484 




PCT/DE03/01002 



Die Position des Massezentrums des Korpers ist in einem 
extemen Koordinatensystem X± gegeben. Ein zweites 
Koordinatennetz xi ist mit dem Korper verbunden, so dass seine 
5 Oberf lachenpunkte im zweiten System gegeben sind. 

Fur eine dextrale Scherung hat das deviatorische externe 
Kraft f eld die Form F ex tr = fext/ wobei 



F = 

K ext 



0 



dx 2 



(Gl.29) 

0 0 

10 (Abb. 08a), Dieses Feld tritt mit einem Einheitskorper (dem 
thermodynamischen System) mit isotropen Materialeigenschaf ten 
in Wechselwirkung (Abb. 08b). Daher setzen in einem Kontinuum 
eine unendliche Zahl von Kraf tvektoren auf einer unendlichen 
Zahl von Punkten auf der Systemoberf lache an. Daher muss man 

15 den durchschnittlichen Ansatzpunkt des durchschnitt lichen 
Kraftvektors f av auf dem Korper finden. So, wie sich das 
Problem in Abb. 08 darstellt, ist f oc sin9, wobei 8 der Winkel 
zwischen r und der Koordinate x x ist. Der Durchschnitt fur 
einen Quadranten ist dann 

20 

psin9(/9 2 

* = -. (G1.30) 

n n 

2 

Die durchschnittliche Kraft ist ein Vektor der Form f av = [2/n, 
0] , der Positionsvektor seines Ansatzpunktes P K bildet mit x x 
den Winkel K = ± 39,54° (Abb. 09) . Wenn der Korper im Raum 
25 festgehalten wird und einem Kraf tf eld wie in Gl.29 ausgesetzt 
ist, hat die Normalkraft die Form 

f n =(fn)n = cosG sinG (G1.31) 
w&hrend die Scherkraft durch 

f s =(f t)t = sin 2 e . (G1.32) 
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gegeben ist. An P K ist letztere daher 

1.405. (G1.33) 



Wenn der Korper auf das angesetzte Kraftfeld reagieren darf, 
ist f 8 (,c) ausgeglichen: da Koharenz zwischen System und Umgebung 
5 besteht, ist ein ausseres Ungleichgewicht unmoglich. Daher kann 
das System nicht frei rotieren, sondern es treten die 
Oberf lachenbindungskraf te mg in Erscheinung, die das Rotieren 
verhindern. Daher wird |f a ool von alien |f s l an alien Punkten 
abgezogen, wodurch das Vorzeichen des effektiven f s in einigen 
10 Gebieten umgekehrt wird. Die effektive Scherkraf t-Grofie an 
alien Punkten P ist daher 

M.-vt.-M. (G1 - 34) 

so dass die gesamte effektive Kraft an P(0) 



f eff =W f n- " (G1.35) 



15 ist 



Zu Beginn gab es zwei Punkte ohne Kraft komponente senkrecht zum 
Radius, an 9 = 0 und 8 = it (Abb. 08b) . Durch die Subtraktion von 
f s(K , entstehen stattdessen zwei neue mit Winkelabstand 9 = ± k 

2 0 auf beiden Seiten der Xi-Koordinate . Da entlang dieser 
Richtungen nur Norma lkomponen ten aktiv sind, werden sie als 
Eigenrichtungen interpretiert (kontrahierend : c, extendierend : 
e) , die nicht mehr orthogonal sind. Sie schliefien miteinander 
die Winkel n - 2k = 100.92° und 2k = 79.08° ein. Jedoch wurden 

25 oben in den Gl.31 und G1.32 die Einheitsvektoren n und t ± n 
ohne weitere Erlauterung benutzt; implizit wurde damit jedoch 
eine spharische Form fur das Volumenelement vorausgeset zt , da 
es zwischen n und r keinen Unterschied gibt. Es ist jedoch 
immer noch nicht moglich, das Drehmoment fur einen spharischen 

30 Korper auszugleichen, da \ 2n fsieffjx^ d0 * 0 . Daher wird 
angenommen, dass die Annahme eines spharischen Korpers 
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inkorrekt ist. Die Winkelbeziehung zwischen den Eigenrichtungen 
legt nahe, dass sowohl das effektive Kraft f eld £ ef f fur einfache 
Scherung als auch die Form des Systems elliptische 
Eigenschaften haben. Daher ist es notwendig, die Form der 
Ellipse zu berechnen, die mit f e ff mit den Eigenrichtungen wie 
oben im Gleichgewicht ist, und die die Bedingungen von Gl.18 
erfullt . 



Form von Kraftfeld und Volumenelement 

Im Beispiel der reinen Scherung (RS) hat der 
Feldeigenschaf tstensor die Eigenschaften tr Prs = 0 und 
det Frs a ± l. G1.2 9 auf einfache Scherung (ES) angewandt 
weicht davon ab. Jedoch ist der Germ in G1.33 nicht Teil von 
G1.29. Daher wird angenommen / eiavtf: die Feldmatrix fur einfache 
Scherung nach der Siriatr^^v*: v - Eigenschaf t det F ES = *± u l 
hat, was anzeigt, dass alle Pi&^se, die vom deviatorischen Feld 
represent iert werden und ins System eintreten, es auch wieder 
verlassen. Daher werden das Kraftfeld fur die Bedingungen der 
elastischen einfachen Scherung und die Forme igenschaf ten des 
Systems hier modelliert, indem ein orthogonales Feld der Form 

Mo -°J 

angenommen wird, das im Gleichgewicht mit einem kugelf ormigen 
System ist. Die Hauptscherrichtungen sollten jedoch mit den 
Punkten unter den Koordinaten zusammenf alien, weshalb das Feld 
urn 45° sinistral gedreht wird. Sowohl Kraftfeld als auch 
Radiusfeld werden dann mit Hilfe einer Trans format ionsmatrix T 
umgewandelt, so dass sie elliptische Eigenschaften erhalten, 
wobei p die Positionsvektoren der Oberf lachenpunkte im nicht - 
transformierten Zustand sind: 



[cot* 0 1 
"[ 0 tanicj' 



(G1.37) 
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15 



r o cotKi r o r„i 

Tp = r eD =[r n cos6 r^sinG], (G1.39) 

und 

FdcvP = fdev=fciSine T^cosG]. (Gl.40) 

5 

Die Eigenrichtungen des elliptischen Radiusfeldes r e ii sind 
daher orthogonal und identisch mit den Hauptachsen der 
resultierenden Ellipse, die Querschnittsf lache A/n = T 1X T 2 2 = 1 
hat, wahrend die Eigenrichtungen v e und v c des Kraftfeldes wie 
10 gewunscht nicht -orthogonal sind (Abb. 10), 



/c * c [±tanj 



(G1.41) 



Kinematic . : "reaction Scherung 



Es trifft sich gunstig, dass alle die Gleichgewichtsbedingungen 
betreffenden Gleichungen nur die transf ormierten Gleichungen 
des Beispiels der reinen Scherung sind. Das Kreuzprodukt 
vereinfacht sich zu 

2 0 ^ x r en | dQ = T U T 22 f(cos 2 G - sin 2 G) rfG = 0 , (Gl . 42 ) 

d.h. die elliptischen Eigenschaf ten von Radiusfeld und 
Kraftfeld kurzen einander raus. Das lineare Moment ist 

ffdev ' r eii d ® = ( T u + T 2i) fcosGsinG dQ = 0 . (Gl . 43 ) 

T? t + Tf 2 >2, ist die maximale GroSe des Skalarprodukts entlang 

25 der Eigenrichtungen >0. Dies wird als Artefakt der elliptischen 
Form angesehen, da die Radii entlang der Eigenrichtungen wegen 
der Transformation durch T nicht Einheitslange haben. Wenn 
dagegen das Skalarprodukt des inversen Positionsvektors r e ii" 
1 = T _1 r durchgefuhrt wird, 

3 0 ^v)! dB = jr* ' fdcv dQ = j 2cosGsinG dQ = 0 (Gl . 44 ) 
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ist die Gleichgewichtsbedingung eingehalten. Die elliptischen 
Eigenschaf ten von Kraftfeld und Radiusfeld heben einander auf , 
was den Bezug einer normalen deviatorischen Kraft f(n)dev pro 
Einheitsradius zur Folge hat, deren maximale GroSe 1 ist. Gl.44 
zeigt die normalisierte Beziehung von | f n <dev> | zu |r el i|; diese 
Normal isierung von f n mit Bezug zu r ist notwendig (vgl . 
Gl.23). fop (vgl. Gl.19-Gl.23) ist daher auch ein 
Einheitsvektor. Gl.42 und Gl.44 zusammen halten G1.18 ein. 
Einige Vektorgrofien sind in Abb. 11 gezeigt. 



Wenn t JL n, ist t kein Tangent envekt or fur die Ellipse, 
sondern definiert als ein Einheitsvektor senkrecht zum Radius 
und unabhangig von der Orient ierung der Ellipsenoberf lache . 

Die GroSe der Scherkraf tkomponente ist durch 

■«-*■„ sin 2 e-r 22 cos 2 e (01.45) 

gegeben. Wenn Gl.53 fur separate Sektoren (z.B. von der 
kontrahierenden Eigenrichtung in eine von beiden Seiten zur 
nachsten extendierenden Eigenrichtung, oder von 0 bis k und von 
k bis V2) integriert wird, zeigt sich, dass es ein 
Ungleichgewicht zwischen dextralen und sinistralen Scherkraften 
auf der Systemoberf lache gibt . Jedoch ist G1.45 keine 
Gleichgewichtsbedingung, und die Scherkraf te alleine mussen 
nicht ausgeglichen sein. Da die Krafte fur die Ellipse 
ausgeglichen sind (Gl.42), bedeutet ein Ungleichgewicht fur 
fs(dev) keine freie Drehung. Das Ungleichgewicht ist daher 
dahingehend interpretiert , dass das Kraftfeld permanent 
reorientiert wird durch eine Transformation der inneren 
Koordinaten x± gegen die ausseren X± urn einen Winkel X, 

VWo +f S (d« 0 »o=tanX; X = 28.83<» (Gl.46) 

wobei die VektorgrdSen unter der Wurzel fur die integrierten 
Sunvmen iiber benachbarte Sektoren stehen. Die Wurzel wird 
gezogen, weil alle Winkelterme in Gl.45 Quadrat -Terme sind. Da 
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die Eigenrichtungen real sind, ist der Winkel X stabil . 

Es ist zweckmaSig, die linke Seite in Gl.46 als Tangens-Term zu 
interpretieren, weil es keinen Grund gibt, warum der Term auf 
5 Werte < 1 beschrankt sein soil- Wenn LHS -> «>, X -> 90° , was 
die Bedingung fur externe Rotation ware. Im selben Moment 
sollte die charakteristische Gleichung degenerieren, d.h. nur 
eine Losung haben, worauf die Eigenrichtungen imaginar werden.) 

10 Die Erfindung beinhaltet, dass eine Scherkraft einen Punkt P 
nach P 1 verschiebt, weg von c und hin zu e. Eine Scherkraft 
wird immer eine lokale Dehnung proportional zu ihrer GroSe 
bewirken. Dieser Effekt ist additiv von Punkt zu Punkt. 
Derselbe Effekt gilt fur elastische einfache Scherung, ausser 

15 dass die Eigenrichtungen nicht mehr orthogonal sind. Das 
Ergebnis ist jedoch ahnlich zu, dem von '61.2 8 (di* integriert 
uber je einen Quadranten 1 ergab) , nur ist die ■ • : A S.c-h s . uber 
die integriert £ wird, groSer als ein -i^^^chendea- 
Kugelsegment . Das Ergebnis pro Quadrant ist daher im Vergleich 

20 zur Kugel (urn den Faktor) T n T 22 /2 = 1,018. Wie im Beispiel der 
reinen Scherung erklart (nach G1.28), zahlt dieses Ergebnis 
doppelt, also w B = 2,036. 

w Q ist der dilatierende Effekt der Scherkrafte auf den Radius - 
25 r e parallel zu e. Der Kurzungsf aktor fur den Radius r c entlang 
von c wird durch Gl.22 gegeben. Wenn v c und v e normalisierte 
Eigenvektoren sind, konnen sie als die Radien eines 
Einheitskreises aufgefasst werden, der dem monoklinen Feld Pr = 
f ausgesetzt ist, und das Verschiebungsf eld kann ausgerechnet 
30 werden. Ein numerisches Beispiel unter Verwendung von Gl.22: 
r 0 = 1, fo = 1/ und das Volumen V Q /n = 1. Nach Belastung auf die 
operative Feldstarke -1/5 sind die GroSen entlang von c 
Af = f op + fn(dev) = -2/5, r c = |v c l = 0.670. Entlang von e sind 
sie f 0 p + fn(dev) = -1/5 + 1/5 = 0, r e = 1.000; aber 
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Wa = 2.036/5 = 0.407, daher r e = IvJ = 1.503. Wenn S der 
Verschiebungsf eld-Eigenschaftstensor ist , 

|vJvJ = detS = 1.0155 (Gl.47) 
was besagt, dass in diesem Beispiel das Volumen eine reale 
5 Expansion von 0,75% erf ahrt . 

Geometrische Eigenachaften einer Deformation nach einfacher 
Scherxmg 

10 

Die Koordinaten Xi werden gegenuber X± urn 1 in der 
Scherrichtung rotiert. Im ausseren Rahmen X± ergibt sich 
folgendes Bild (Abb. 12) : die extendierende Eigenrichtung e 
liegt bei 10,71°, die kontrahierende Eigenrichtung c bei 

15 111,63°. Die Winkelhalbierenden der von c und e 
eingeschlossenen Winkel sind. parVliiu r.u den x±- Koordinaten und 
. . markieren die Richtunge^ s v^^aler Scherung : Der 
Orient ierungsvektor, der den grSSeren Sektor bei 61,18° 
halbiert, markiert eine Scherrichtung senkrecht zu sich selbst, 

20 die R-Flache liegt daher bei -28,83°; der Orient ierungsvektor , 
der den kleinen Sektor bei 28,83° teilt, markiert die P- 
Scherflache bei 61,18° (Abb. 12). 

Ein wichtiger Bestandteil der Erfindung ist die Ermittlung des 
2 5 Kraftfeldes, das die Deformation verursacht, und des 
Verschiebungsf eldes. Jedoch muss das Verschiebungsf eld 
unabhangig von der Def ormationsart - elastisch oder plastisch - 
die Eigenschaf ten des Kraftfeldes widerspiegeln, selbst wenn 
die dabei entstehenden Phanomene sich ganz erheblich in ihrer 
30 Natur unterscheiden. Daher ist es gerechtf ertigt , die 
Eigenschaf ten des berechneten Kraftfeldes mit 

Gefugeeigenschaf ten zu vergleichen, die in plastisch verformten 
Gesteinen zu beobachten sind. X 1 ist die Scherzonengrenze 
(Abb. 12a) . Von den Eigenrichtungen wird e mit der S-Flache in 
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SC-Gefuge korreliert (Berthe et al . 1979) und mit der Richtung, 
parallel zu der die Hauptanisotropie der Kristalle sich 
ausrichtet. Der vorausgesagte Winkel von ca.ll° ist in guter 
Ubereinstimmung mit der Schiefe von Gef ugediagrammen fur 
monomineralische Scherzonen mit Mineralen mit nur einer gut 
ausgebildeten Scherrichtung (Eis, Glimmer, Olivin) . 

Die kontrahierende Eigenrichtung c ist in Hochtemperatur- 
Tektoniten nicht gut entwickelt. Die Richtung kann durch Fehlen 
von Druckschatten oder durch minimale Glimmer- Ausrichtung 
entlang der Oberflache von Feldspat-Porphyroklasten erkannt 
werden. Wenn die Mylonite jedoch herausgehoben werden, 
entwickeln sie haufige Klufte, die den Lagenbau mit ca. 70-80° 
gegen die Scherrichtung geneigt schneiden. Sie sind offenbar 
von der elastischen Energie verursacht, die wahrend der 
plastischen Verfcrmung im Gestein gespeichert wurde, und die 
freigesetzt wixi/ wenn der UmschlieSungsdruck nicht r ; mehr 
ausreicht, 6* : &%-<3kife\£&in zusammenzuhalten. Das hier vorgestellte 
Verfahren erlaubt den Schluss, dass Risse sich parallel zu c 
of f nen. 

Die Scherflache bei 0 = -28,83° ist die Riedelsche Flache R 
(Riedel 1929) oder die C-Flache in SC-Gefuge (Berthe et al . 
1979) . Die P-Flache ist bei naturlicher plastischer Verformung 
meistens unterdruckt, wurde aber in Experimenten beobachtet 
(z.B. Merzer & Freund 1975), bildet auch manchmal eine schlecht 
ausgebildete Scherrichtung (Abb. 12b). Sie tritt starker in 
Erscheinung, wenn das Materialverhalten nahe der plastisch- 
bruchhaften Grenze lag. Die R-Flache (synthetisch, hier 
dextral) und die P-Flache (antithetisch, hier sinistral) sind 
nicht mechanisch gleichwertig. Die R-Flache ist wohl die erste, 
die zu flieSen anfangt, wenn plastisches Verhalten einsetzt. 
Abb. 13 zeigt das FlieSfeld. Der Schersinn auf der R-Flache ist 
synthetisch zur Gesamt scherrichtung, aber die Flache selbst 
rotiert gleichzeitig antithetisch nach e; bei diesem Vorgang 
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wirci sie progressiv gestreckt, bis sie zerfallt. R- und P- 
Flachen sind Dehnungsgleitf lachen (stretch faults) . 

Energetik der Deformation Jbei elastischer und plastischer, 
reiner und einfacher Scherung 

Elastische Verformung 

Der Zustand, in dem das System sich befindet, wird vom 
operativen Kraft f eld bestimmt, das als MaS fur die Feldstarke 
dient. Es kann genutzt werden, die Def ormationstypen 
miteinander zu vergleichen. Im elastischen Bereich sind 
operatives und deviatorisches Feld wahrend der progressiven 
Belastung proportional zueinander; sie benotigen daher beide 
Arbeit, die das elastische Potential reprasentiert . Da die 
.iumme von f n <dev) von 0 nach n Null ist, ist die von 
Normalkraften geleistete Arbeit nur die von fop/ 



fur r 0 = f o = Af op = 1 (vgl . Gl.22) . Im isotropen Fall ware dies 
der einzige Term, der in 2D zu berucksichtigen ware; in 
absoluten Zahlen ist dies von alien Moglichkeiten der Fall der 
groSten Arbeit. Diese Arbeit aussert sich in isotroper 
Kontraktion. Von Scherkraften geleistete Arbeit verursacht eine 
gerichtete Dilation. Gl.28 integriert von 0 nach p/2 fur einen 
Quadranten, mal 4 fur alle Quadranten, ergibt 8,000; die 
Gesamtarbeit fur elastische reine Scherung betragt daher 
Wps(el) = 8 - 271 = 1,717. 

In elastischer einfacher Scherung ist die GroSe von f op 
dieselbe wie oben. Die Scherarbeit (G1.47) muss fur einen 
Sektor von c nach e integriert und mit 4 multipliziert werden, 
urn die Summe liber alle Sektoren zu erhalten. Die Gesamtarbeit 
bei elastischer einfacher Scherung ist daher ws S < e i) = ~ 271 + 




(G1.48) 
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8,144 = 1,861. Bei gleicher operativer Feldstarke erfordert 
elastische einfache Scherung daher 8,4% mehr Energie als reine 
Scherung. 

5 

Plastische Verformung 

Der Mechanismus, durch den eine plastische Verformung bewirkt 
wird, ist nicht Gegenstand der Erfindung. Es wird angenommen, 

10 dass das Material kontinuierlich bleibt, und dass ein homogener 
Zustand des Feldes f f ortexistiert . Der Ubergang von 
elastischer zu plastischer Verformung (der FlieSpunkt) wird 
durch einen Wechsel von reversiblem zu irreversiblem Verhalten 
charakterisiert . Daher erreicht das operative Feld eine stabile 

15 Feldstarke und erfordert keine weitere Arbeit. 

Die gesamte weitere Arbeit, die plastische Verf ormvu:' ^ : .^rg^e;, > ' 
wird daher nur durch Scherkrafte geleistet. Da 'die- ^rr-..ifeien * 
verursachte Dilation ein elastisches Phanomen ist, dessen 

20 Aufbau am FlieSpunkt aufhort, wird die Energie nicht im Volumen 
des Systems gespeichert, sondern entlang seiner Oberflache 
dissipiert. Da das System sich daher passiv verhalt, gilt das 
Prinzip der doppelten Arbeit (vgl . Diskussion nach Gl.28) nicht 
mehr; also ist die benotigte Scherenergie nur halb so groS wie 

25 fur eine elastische Verformung. Dissipierte Energie ist 
verlorene Energie und muss nicht ausgeglichen sein; da aber 
eine freie Drehung des gesamten Systems in einem Feld mit 
realen Eigenrichtungen nicht moglich ist, - die die nach wie 
vor existieren wegen der f ortgesetzten Existenz des elastisch 

30 gespannten Zustands - muss das Drehgleichgewicht nach wie vor 
ausgeglichen sein. 

Fur einfache Scherung ist das Feld f orthorhombisch, und 
sinistrale und dextrale Scherung gleichen einander aus, also 
35 ist eine externe Rotation des Systems nicht moglich. Bei der 
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Betrachtung der Arbeit sind die sinistrale und dextrale 
Scherarbeit zu summieren, daher ist die bei plastischer reiner 
Scherung geleistete Arbeit halb so groS wie bei elastischer 
Belastung, also wp S < p i) = 4. 

Bei einfacher Scherung sind die Verhaltnisse etwas komplexer. 
G1.45 von c nach e integriert fur beide Sektoren ergibt 
verschiedene Werte : 1,229 dextrale Scherarbeit in den groSen 
Sektoren, und 0,733 sinistrale Arbeit in den kleinen Sektoren. 
Das Ungleichgewicht fur das gesamte System betragt 1,202; die 
gesamte auf das System geleistete Scherarbeit betragt 4,144. Da 
die grofien Sektoren das Ubergewicht haben, lost das 
Ungleichgewicht einen plastischen dextralen Gesamtfluss aus. 
Die Auslosung eines dextralen Ungleichgewichts fuhrt zu einem 
Ungleichgewicht des Drehmoments fur das gesamte System. Da die 
Form des Systems elliptisch ist und dis J.ange Achse des Systems 
im kleinen Sektor mit sinistra',? vvferung ist, ist --das 
Drehmoment des . gesamten Systems atonal zur GroSe >^%ber 

umgekehrt zum Drehsinn der plastischen Scherung, also 
sinistral. Die dextrale plastische Scherung kostet Energie, die 
dissipiert wird, d.h. es ist eine innere dextrale Scherung. Die 
Rotation des Systems aufgrund des Drehmoments ist eine aussere 
sinistrale Scherung, die keine Energie kostet. Da innere und 
aussere Rotation einander ausgleichen, sind die Eigenrichtungen 
immer noch real. Daher ist also die dissipierte Energie bei 
einer plastischen dextralen einfachen Scherung durch eine 
externe sinistrale Rotation minimiert, und die so gesparte 
Energie ist proportional zum Ungleichgewicht; also ist w es . 
pX = 4,144 - 1,202 = 2,933. 

Bei gleicher operativer Feldstarke kostet eine plastische 
einfache Scherung daher 26,7% weniger Energie als reine 
Scherung. Energetische Unterschiede dieser Grofienordnung sind 
von Experimenten bekannt; sie treten systematisch auf als 
Funktion des Verschiebungsf eldes und sind materialunabhangig 
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(Quarz, Salz, Gips; Franssen & Spiers 1990; Franssen, pers. 
Mitt . ) . 

Das hier vorgestellte Verfahren ist vollstandig mit der 
5 Potentialtheorie und thermodynamischen Prinzipien vereinbar. 
Die Erfindung bietet eine einfache, unmittelbare Korrelation 
von Kraftfeld zu Verschiebungsf eld an, die beide Vektorfelder 
sind und von Tensoren kontrolliert werden; aus mathematischer 
Sicht ist das systematischer als das bisherige Verfahren. Das 
10 neue Verfahren macht ausserdem Vorhersagen moglich, die das 
bisherige Verfahren nicht gestattet, weil es z.B. aus einem 
gegebenen Spannungszustand nicht die Ableitung eines bestimmten 
Verschiebungsf eldes erlaubt . 

15 Ein viskoser Def ormationsschritt fur ein Gas oder eine 
Flussigkeit kann ssrlegt werden in einen zeitunabhangigeri, 
reversiblen, ^.■.^:;p:ben Belastungsschritt .und einen 
zeitabhangigen, ^reversiblen, dif fusionskontrollierten 

Relaxationsschritt . Das neue Verfahren ist auch geeignet zur 

20 Modellierung des Verhaltens von viskosen Flussigkeiten, sofern 
Zeitraume betrachtet werden, die im Vergleich zur Relaxations - 
geschwindigkeit kurz sind. 

25 Anwendungsbei spiel fur eine elastische Deformation eines 
Werks tucks mit gegebener Form unter gegebenen Grenzbedingungen 

Das vorgestellte Verfahren basiert auf Potentialberechnungen. 
Das thermodynamische Kontinuum ist dabei kein Kontinuum von 

3 0 Punkten, sondern ein Kontinuum von Systemen. Ein beliebiger 
Punkt Q im mit verteilter Masse gefullten euklidischen Raum 
dient als Sitz eines thermodynamischen Systems; seine 
physikalischen Eigenschaf ten - Dichte, chemische 

Zusammensetzung, Temperatur, Materialeigenschaf ten, 

35 Orientierung der Anisotropie - und zusatzlich die der externen 
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Grenzbedingungen sind dann Ortsf unktionen von Q. 
MaSstabunabhangigkeit stellt dann sicher, dass die extensiven 
GroSen pro Einheitsmasse gultig sind. Zwei unendlich nah 
beieinander liegende Punkte Qi und Q 2 represent ier en zwei 
unterschiedliche thermodynamische Systeme Vi und V 2 mit einer 
finiten Grofie, die zum groSten Teil uberlappen. Ihre 
physikalischen Eigenschaf ten beziehungsweise die 

Grenzbedingungen konnen Gradienten im euklidischen Raum 
unterworfen sein. Die Systeme konnen so definiert sein, dass 
sie Einheitsmasse, Einheitsvolumen oder Einheitsradius haben, 
die allesamt finit sind. Im Unterschied zum bisherigen 
Verfahren bietet die hier vorgestellte Erfindung eine seiner 
inneren Logik folgende finite Strecke - den Radius des 
thermodynamischen Systems - als Basisstrecke fur eine 
Wellengleichung an. Die raumlichen Eigenschaf ten der 
phv^ikalischen GroSen konnen daher durch einen Ansatz unter 
;t.^;?V.dung von Fourierschen Serien als Funktion der 
^^t^bedingungen modelliert werden. 

Das gesamte Kraftfeld ist auf gegliedert in das isotrope 
Kraftfeld f op und das deviatorische Kraftfeld f de v (vgl . Teil I, 
Koenemann 2000) . Letzteres kann weiter auf gegliedert werden in 
die Normalkomponente f n <dev> und die Scherkomponente f 9 <dev) . 
Normalkrafte werden den ihnen zugeordneten Radiusvektor je nach 
Vorzeichen verkurzen oder strecken. Scherkrafte dagegen haben 
immer dehnende Wirkung, was den Effekt einer zusatzlichen 
Streckungskomponente parallel zur extendi erenden Eigenrichtung 
e hat, die parallel zu x x ist. Nach Gl.28 (Koenemann 2000) 
wurde f estgestellt , dass die von Scherkraften geleistete Arbeit 
doppelt gilt, weil der eine Teil der Streckung in x x vom System 
selbst, der andere jedoch von der Umgebung geleistet wird. Der 
vom System geleistete Teil ist eine Funktion der 
Belastungskonf iguration (BelastungsgrdSe, externe 

Grenzbedingungen, z.B. reine Scherung) . Die von der Umgebung 
verursachte Scherdilatanz gibt es in einem perfekten Kontinuum 
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durchaus, aber ihre Grofie muss auf Grenzf lachen gegen den 
freien Raum Null erreichen. Wenn alle anderen Komponenten 
konstant gehalten werden, ist es diese Komponente, die als 
Funktion des Ortes innerhalb eines Korpers variieren kann und 
5 damit eine Funktion der Form des Korpers ist. 

Das 2D-Modell besteht aus einem rechteckigen Festkorper. Es 
wird oben und unten mit einer konstanten Kraft belastet, 
wahrend die Seitenf lachen frei sind. Die Last wirkt vertikal 

10 parallel zu y {df y /dx = 0) . Die belasteten Flachen (y = 0 , y = 
1) sind wegen der Reibung an den Stempeln in ihrer Freiheit 
behindert und konnen sich in x nicht ausdehnen; an alien 
anderen y sind sie frei, Daher ist die Belastung im gesamten 
Korper dieselbe, es gibt keinen Gradienten, einschlieSlich der 

15 Punkte auf den Seitenf lachen. Daher sind die GroSen des 
operativen Feldes f op und der Normal en -Komponente ' £n(^- ■ i*n 
gesamten Korper monoton. Sie werden hier nicht weitfer iw-^r-'ca^t- 
betrachtet . 

2 0 Die Grenzbedingungen fur die dilatierende Scherung T variieren 
im Korper auf zweierlei Weise: (a) Die Bedingung, dass die 
belasteten Flachen sich in x nicht ausdehnen durfen, erreicht 
ihre maximale Wirkung auf den Flachen y = 0 und y = 1, ist aber 
gegen die Mitte des Korpers hin zunehmend unwirksam, was eine 

2 5 progressive Ausdehnung in x gestattet (hier genannt T x ) ; und 

(b) der Kontinuum-Scheref f ekt (T 2 ) erreicht ein Minimum - auf 
den freien Seitenf lachen, weil es keine Umgebung mehr gibt; auf 
den belasteten Flachen oben und unten, weil Ausdehnung durch 
Bedingung (a) verhindert wird - sie steigt aber gegen die Mitte 
30 des Korpers an. Die Richtung der abnehmenden Strenge fur die 
Grenzbedingungen ist gegen die Mitte fur T x und nach alien 
Seiten hin fur T 2 . Also haben die beiden Grenzbedingungen 
entgegengesetzte Wirkung und Vorzeichen (Abb. 14, Pfeile im 
vierten Ring fur die Bedingungen 1 und 2) . Sie mussen daher 

3 5 separat modelliert werden. 
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Im statisch belasteten Zustand leisten System und Umgebung 
miteinander Arbeit, die nicht gemessen werden kann, weil in der 
Thermodynamik nur die mit Zustandsanderungen verbundene Arbeit 
gemessen werden kann. Zeitabhangige Prozesse werden nicht 
beobachtet, System und Umgebung sind im Gleichgewicht 
miteinander, und alle Flusse sind ausgeglichen. Die zur 
Diskussion stehende Grofie T = 3i + T 2 ist die von Scherkraften 
geleistete Arbeit, eine skalare Grofie, die eine Volumenanderung 
bewirkt, die ebenfalls eine skalare GroSe ist. Es ist daher 
moglich, das Modell mit Hilfe einer Laplace-Gleichung 
anzugehen . 

Ti wird zuerst behandelt. Die Verteilung einer skalaren Grofie 
ti in der Region x(0->b) y(0->d) mit einem Maximalwert ti - -M fur 
alle Punkte y = 0 und ti = 0 fur alle Pt.mfcte x . - 0, x = b .und 
y = d ist 



(Boas 1983:546, ihre Gleichung 2.17), wobei M ein die absolute 
Grofie bestimmender Faktor, Jb die Weite des Korpers in x, d 
seine Lange in y, und m = 2 n - 1 ist (jn ist immer ungerade) . 
Der umgekehrte Fall ist der, bei dem t 2 = M fur alle Punkte y = 
d, und M = 0 ist fur alle Punkte x = 0, x = b, und y = 0, 




mnsinh 



(G1.49) 



b 




mnsinh 



(G1.50) 



Der gemeinsame Effekt ist T x = ti + t 2 . Entlang x = b/2 bildet 
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T x eine nach oben gewolbte Kurve mit einem Minimum Ti = -M fur 
y = 0 und y = d und einem Maximum (M<0) fur y = d/2 . Entlang y 
= d/2 bildet Ti eine nach unten gewolbte Kurve mit einem 
Maximum Ti fur x = 0 und x = Jb und einem Minimum bei y = d/2 
5 (welches der maximale Punkt auf x = Jb/2 ist) . 

Fur T 2 sind die Grenzbedingungen wie f olgt : M = 0 fur alle 
Punkte x=0, x = Jb, y=0 und y = d, und der Graph steigt nach 
innen an. Dies geschieht durch 

10 

T 2 =bsin^dsin?2- (G1.51) 

Die Summe T = T x + T 2 wird als das Potential fur Dehnungsrisse 
interpretiert (in Abb. 15-Abb. 18 jeweils oben). Es ist von 
15 besonderem Interesse, berauszuf inden, wo die Gradienten fur 
Gl . 3 am groSten siu-S.. -L:r?.\v^r 

*± = -AM±^ b - A b ) U } (G1.52) 



20 



n=1 bsinh- 

b 



und 



^L = 7CC05 ™ rfjfcS; (G1.53) 
ax ft d 



die Grofie des Gradienten wird durch 



. mux 
\mnsin- 



^ = 4A/f k * t*-! b - (Gl . 54 ) 



25 und 



7i sin — dsxn — 
o 1 2 ft d 



dx 2 



L = 2 (G1.55) 
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gegeben. Die Ableitung d^T/dx? (die Summe von G1.546 und G1.55) 
der Oberflache T(b, d) wird als das Potential fur Scherrisse 
interpretiert (Abb. 15 -Abb. 18 jeweils unten) . Die Ableitung 
5 d^T/dy* liefert dieselben Ergebnisse in dem Sinn, dass die 
Morphologien der beiden Ableitungen identisch sind, jedoch sind 
ihre absoluten GroSen von den gewahlten absoluten GroSen von Jb 
bzw. d abhangig. Fur realistische Ergebnisse muss das Modell 
mit Hilfe von beobachteten Werten geeicht werden. 

10 

Der Ansatz liefert alle Eigenschaf ten der Spannung als Funktion 
des Ortes innerhalb der Gultigkeitsregion Jbd. Auf ahnliche 
Weise sollte es moglich sein, die Verschiebungen aus der 
Arbeit sglei chung herzuleiten und die Form des deformierten 
15 Korpers zu berechnen. 

* ■ ^ 

20 Abb.l5-Abb. 18 zeigen die oberen Teile das Potential fur 
Dehnungsrisse T, die unteren Teile zeigen das Potential fur 
Scherrisse SPTdx*. Die Breite von 10 Entf ernungseinheiten ist 
dieselbe in alien Abbildungen. Die Farben deuten relative 
GroSen an (lila = Minimum, rot = Maximum) . 

25 

Fur dicke, kurze Saulen (Abb. 15, Hohe 2 Einheiten) ist das 
Zentrum dynamisch tot, da es unter starker Kompression steht . 
Das Dilationspotential ist ausschlieSlich entlang der freien 
Oberflachen zu finden, was nahe legt, dass ein bis an die 
30 elastische Grenze belastetes sprodes Material oberf lachlich 
absplittern wird. Das Scherrisspotential ist weiter innerhalb 
des Korpers, was die Oberflache weiter destabilisiert , es ist 
aber sehr niedrig im Zentrum. Die Spannungsverteilung wird 
stark von der Bedingung kontrolliert , dass die belasteten 
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Oberflachen oben und unten nicht entlang der Stempel gleiten 
durfen. 

Wenn die Hohe des Korpers nur wenig groSer gewahlt wird 
(Abb. 16, Hohe 5 Einheiten) , befindet sich das Dehnungsriss- 
Potential immer noch nahe an der Oberflache, aber das 
Scherrisspotential ist sehr schnell ins Werkstuck-Innere 
gewandert. Die beiden Maxima fallen nicht mehr aufeinander. 

In einem wurf elf ormigen Korper (Abb. 17, Hohe 10 Einheiten) 
fallen das Dehnungsriss- und das Scherrissmaximum wieder 
aufeinander, aber diesmal im Zentrum des Werkstucks, 
dahingehend interpretiert wird, dass jegliche Form des 
Materialnachgebens hier anfangt, sei es bruchhaft oder 
plastisch. 

Schliefilich wird ein langer, schmaler Korper betracx^BU 
(Abb. 18, Hohe 30 Einheiten). Das Dehnungsrisspotential mMfe 
im Zentrum des Werkstucks, aber es breitet sich entlang der 
langen Achse aus und nimmt nur kurz vor den belasteten Flachen 
ab. Das Scherrisspotential ist von der Mitte in Richtung auf 
die belasteten Flachen hin gewandert und bildet zwei Maxima, 
die gegen die Kanten hin auslaufen, an denen die belasteten und 
die freien Flachen aneinander stoSen und eine Grenzbedingungs- 
Diskontinuitat existiert. Die beiden Auslaufer zu den Ecken und 
die belastete Flache umschlieSen eine dreieckige Region, die 
ausserdem nur geringer Dehnung ausgesetzt sind; diese Dreiecke 
sind praktisch tot. Die Berechnungen erlauben die 
Interpretation, dass bruchhaftes Verhalten aufgrund von 
Dehnungsrissen im Zentrum des Werkstucks beginnen sollte; die 
Risse sollten sich entlang der langen Achse des Versuchskorpers 
fortsetzen (was letztendlich die lokalen Grenzbedingungen 
innerhalb des Korpers in gewissem Grade andert) , bis der 
Dehnungsef fekt schwindet; aber genau dort erreicht das 
Scherrisspotential sein Maximum, weshalb zu erwarten ist, dass 
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der Bruchmodus von Dehnungsriss auf Scherriss ubergehen und die 
Risse sich gegen die Kanten des Korpers fortsetzen sollten. Das 
vorhergesagte Materialverhalten wird in der Tat in sproden 
Materialen beobachtet, z.B. in Keramik (Winter 1992). 



Schlus s £o lgerung 



zur 



Es wurde gezeigt, dass das neue Verfahren 
Deformationsermittlung so moduliert werden kann, urn die 
Spannung in einem diskreten Korper darzustellen, der einer 
bestimmten Belastungskonf iguration ausgesetzt ist . Im Vergleich 
zum Aufwand fur das bisherige Verfahren sind die 
Voraussetzungen an technischer Ausrustung und mathematischem 
Gerust urn die hier vorgestellten Ergebra:£3£ zu .erzielen 
ausserst gering. Die beinhalteten Rechnurg-: k \rr^n namlioh auf 
beliebigen Computern durchgefuhrt werden. 

Der Begriff "Computer" ist in keiner Weise einschrankend zu 
verstehen. Es kann sich hierbei urn eine beliebige zur 
Durchfuhrung von Berechnungen geeignete Einheit handeln, 
beispielsweise eine Workstation, einen Personalcomputer , einen 
Mikrocomputer oder eine zur Durchfuhrung von Berechnungen 
geeignete Schaltung. 



Ein wichtiger Vorteil der Erfindung ist, dass kein Punktnetz 
notwendig ist, aber eine Losung kann ohne weiteres fur jeden 
beliebigen Punkt innerhalb des Giiltigkeitsbereiches erlangt 
werden, wobei mathematische Methoden zur Anwendung kommen, fur 
die Losungen von Problemen aus anderen Feldtheorien seit langem 
existieren. 

Insbesondere ermoglicht es die Erfindung, Verformungen ohne 
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Einsatz eines Netzes von Knotenpunkten, wie es beispielhaft in 
Fig. 01 dargestellt ist, vorzunehmen. 

Die Erfindung ermoglicht so eine wesentlich einfachere 
5 Ermittlung von Verformungen. 

Literatur 

Berths , D . , Choukroune , P . , Jegouzo , P . (1979) Orthognei ss , 
10 mylonite and non-coaxial deformation of granites: the 

example of the South Armor ican Shear Zone. Journal of 
Structural Geology 1, 31-41 

Boas, M.L. (1983) Mathematical methods in the physical 
sciences. J. Wiley & Sons, 2nd ed. , 793pp 

15 Bridgman, P.W. (1958) The physics of high pressure. Bell & Sons 
Ltd, London; 445pp. 

Franssen, R.C.M.W., Sr^:^ , . C.J.:- -(1990) Deformation of 
polycrystalline salt in compression and in shear at 250- 
350°C. In: Knipe, R. (ed.) Deformation mechanisms, rheology 
20 and tectonics. Geological Society of London Special Paper 

54, 201-213 

Jeanloz, R., Knittle, E. (1986) Seduction of mantle and core 
properties to a standard state by adiabatic decompression. 
In: Saxena, S. (ed.) Chemistry and physics of terrestrial 
25 planets. Springer Verlag, 275-309. 

Kellogg, O.D. (1929) Foundations of potential theory. Springer 
Verlag, 384pp. 

Malvern, L.E. (1969) Introduction to the mechanics of a 
continuous medium. Prentice-Hall, 713pp. 

30 Merzer, A.M., Freund, R. (1975) Buckling of strike slip faults 
- in model and in nature. Geophysical Journal of the Royal 
Astronomical Society 43, 517-530 



WO 03/081484 




PCT/DE03/01002 



49 



Riedel, W. (1929) Zur Mechanik geologischer Brucherscheinungen. 
Centralblatt fur Mineralogie, Abt.B, 1929, 354-368 

Weber, K. , Juckenack, C. (1990) The structure of the Spessart 
Mtns. crystalline basement and its position in the frame of 
the Mid-European Variscides. In: Mid-German Crystalline 
Rise & Rheinisches Schiefergebirge. Field trip guide to 
IGCP 233 Conference on Paleozoic orogens in Central Europe 
- Geology and Geophysics. Gottingen-Giessen 1990, 101-114 

Winter, W. (1992) Experimen telle Bestimmung der Zugfestigkeit 
sproder Werkstoffe (Glas, Keramik) im Scheiben-Druck- 
Versuch. Diss RWTH Aachen 1992. 



WO 03/081484 




PCT/DE03/01002 



50 



Patentanspruche : 

1. Verfahren zur Ermittlung von Verf ormungen eines Korpers, 
wobei untersucht wird, wie die Einwirkung von Kraften den 
Korper verformt, dadurch gekennzeichnet, dass fur 
einen Teil des Korpers oder seine Gesamtheit ermittelt 
wird, wie in Abhangigkeit von Materialeigenschaf ten und 
ausseren Grenzbedingungen ein Kraf tegleichgewicht zwischen 
irmeren und ausseren Kraften zu Stande kommt, und dass die 
Verf ormungen des Korpers so ermittelt werden, dass alle 
einwirkenden Krafte im Wesent lichen folgender Formel 
entsprechen: 

f ex t + XHayst + IHa + f B (ext) + Hl^gyst) = 0. 

2 . Verfahren nach Anspruch 1 , dadurch 
gekennzeichnet, dass die Materialeigenschaf ten al£ 
Funktion eines Ortes Q ermittelt werden. 

3. Verfahren nach einem oder beiden der Anspruche 1 oder 2, d 
adurchgekennzeichnet, dass 
Kraft feldeigenschaf ten unter Berucksichtigung von extern 
kontrollierten Grenzbedingungen ermittelt werden. 

4 . Verfahren nach einem oder mehreren der vorangegangenen 
Anspruche , dadurch gekennzeichnet, dass die 
gesamte von einem ausseren Vektorfeld geleistete Arbeit 
durch Aufteilung der von Normal- und Scherkomponente eines 
Vektors oder Vektorfeldes geleisteten Arbeit entsprechend 
der Gleichung 

r* f = Jlrxff +|r-f| 2 = |r||f | = const 



berechnet wird. 
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5. Verfahren nach einem oder mehreren der vorangegangenen 



mechanische Wirksamkeit des Radius eines termodynamischen 
Systems im Festkorper berucksichtigt wird. 

6, Verfahren nach einem oder mehreren der vorangegangenen 
Anspruche, dadurch gekennzeichnet , dass innerhalb 
eines Zeitraums, der im Vergleich zur materialabhangigen 
Dif fusionsgeschwindigkeit kurz ist, der Radius des 
thermodynamischen Systems als Hebel berucksichtigt wird. 

7. Verfahren nach einem oder mehreren der vorangegangenen 
Anspruche , dadurch gekennzeichnet, dass die von 
der Scherkomponente der Kraftfelder bewirkte normale 
Streckung f s <ext) und m s(8yst ) im Wesentlichen ^prechend . der 



berechnet wird. 

8. Verfahren nach einem oder mehreren der vorangegangenen 
Anspruche, dadurch gekennzeichnet, dass die 



Anspruche , dadurch gekennzeichnet, dass eine 



Forme 1 




Gesamtverschiebung nach einer 

Arbeit sgleichung, insbesondere der Gleichung 



logarithmischen 




clnr 



ermittelt wird. 



9. 



Verfahren nach einem oder mehreren der vorangegangenen 
Anspruche, dadurch gekennzeichnet, dass die 
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Gesamtverschiebung nach einer aus einer Zustandsgleichung 
hergeleiteten Arbeit sgleichung ermittelt wird. 

10. Computer, dadurch gekennzeichnet ; dass er 
wenigstens ein Berechnungsmittel umfasst, das ein Verfahren 
zur Ermittlung von Verf ormungen eines Korpers durchfuhrt, 
wobei untersucht wird, wie eine Einwirkung von Kraften den 
Korper verformt, wobei in wenigstens einem 
Oberf lachenbereich des Korpers ermittelt wird, ob ein 
Gleichgewicht zwischen inneren und ausseren Kraften 
besteht, dass das Berechnungsmittel 

Gleichgewichtsbedingungen uberprufen kann und dass das 
Berechnungsmittel fur den Fall, dass ein derartiges 
Gleichgewicht besteht, Berechnungen zur Durchfuhrung eines 
Verfahrens nach einem oder mehreren der Anspriiche 1 bis 9 
durchf uhrt . 



